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1. Моделирование 


1.1. Основные понятия и определения 

Моделирование является одной из составных частей аппарата 
системотехники [1, 3, 11]. 

Все то, на что направлена человеческая деятельность, 
называется объектом. 

Модель — это система, имеющая некоторые идентичные 
свойства с оригиналом. 

Замещение реального объекта моделью с целью получения 
информации о важнейших свойствах его называется моделированием. 
Получение информации о реальном объекте осуществляется путем 
проведения экспериментов с его моделью на ЭВМ, АВМ, ЭВМ-АВМ 
или непосредственно на реальном объекте. 

Если эксперименты проводятся на ЭВМ, то их называют 
машинными экспериментами, а сам процесс проведения 
экспериментов имитационным моделированием. Эксперименты, 
которые проводятся непосредственно на реальном объекте, 
называются натурными экспериментами. Эксперименты 
называются полунатурными, если объект, над которым проводиться 
эксперимент помещается в моделируемую или искусственную среду, 
например, при создании новых самолетов эксперименты проводят с 
помощью аэродинамических труб; плазма, полученная в лабораторной 
установке, может быть моделью при изучении звездного вещества. 

Процесс моделирования можно представить в виде цепочки: 

«Описательная модель —математическая схема — 
математическая модель». 

Описательная модель формируется на естественном языке. Под 
понятием «математическая схема» понимают переход от словесного 
описания системы к формальному представлению процесса и его 
функционирования в виде некоторой математической модели 
(аналитической или имитационной). Такой подход позволяет 
рассматривать математику не как метод расчета, а как метод 
мышления. 

При моделировании должны быть рассмотрены три принципа — 
концепции модели: 

• рационализм — аналитические методы на базе фундаментальных 
наук (математика, физика и т.п.); 
• эмпиризм (практика) — экспериментально-статистические методы; 
• прагматизм — целенаправленный деловой подход, в бизнесе — 
инженерный, экономический. 
Рассмотрим основные этапы разработки моделей. 
1. Постановка задачи: необходимость решения самой задачи и 
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моделирования системы; выбор метода решения и определение 
подзадач. 

2. Анализ задачи моделирования системы: выбор критериев 
эффективности процесса функционирования системы; выбор 
возможных методов идентификации; предварительный анализ 
алгоритма модели и результатов моделирования. 

3. Определение требований к исходной информации об объекте 
моделирования и организация её сбора. 

4. Выдвижение гипотез и принятие предположений. 

5. Определение параметров и переменных модели. 

6. Установление основного содержания модели: на этом этапе 
выбирается метод построения модели. 

7. Обоснование критериев оценки эффективности системы. 

8. Определение процедур аппроксимации реальных процессов. 

9. Описание концептуальной модели. 

10. Проверка достоверности концептуальной модели. 

11. Составление технической документации. 


1.2 Общая классификация моделей систем 
Создаваемые модели систем классифицируются по следующим 
основным признакам: 
• по времени: 
= динамические модели: непрерывные, которые описываются 
дифференциальными уравнениями; дискретно-непрерывные 
(разностные), описываются разностными уравнениями; 
вероятностные, построенные на событиях — модели теории 
массового обслуживания; 
" дискретные модели — автоматы; 
• по признаку случайности: 
= детерминированные — модели, отражающие процессы, в 
которых явно отсутствуют либо не учитываются всякие 
случайные воздействия; 


= статистические — модели, отражающие процессы, в которых 
учитываются случайные воздействия; 
• по назначению: экономические, технические, военные, 
социальные, биологические, технологические и т.п. 
• по виду обрабатываемой информации: 
" информационные: справочно-информационные; 


информационно - советующие; экспертные; автоматические; 

=" физические модели: натурные (плазма); полунатурные 
(аэродинамические трубы); 

= имитационные модели; 

= интеллектуальные модели; 


= семантические (логические) модели и т.п. 


1.3 Постановка задачи 

Основными целями моделирования объектов является 

1) изучение его внутренних свойств; 

2) изучение поведения объекта при изменении внешних 
условий влияющих на него. 

Изучение внутренних свойств объекта сопровождается 
уточнением модели объекта и проверкой адекватности модели 
объекту. При этом необходим достаточно большой объем априорных 
(доопытных) данных об объекте, которые получают путем 
наблюдений за его поведением, проведением натурных испытаний и 
т.п. 

Изучение поведения объекта при изменении внешних условий 
влияющих на него проводят на базе модели, не требующей изменения 
в процессе моделирования, и, чаще, данная цель моделирования 
ставится для изучения поведения проектируемых объектов: новых 
приборов, схем, когда получены математические выражения, 
описывающие свойства будущего объекта. 

Моделирование возможно и при объединении обеих целей. 

При достижении целей моделирования проводится 
идентификация объекта моделирования. Существует два типа 
идентификации. 

Структурная идентификация — определение внутренней 
структуры объекта моделирования. Структурная идентификация 
проводится на базе теории системного анализа. 

Параметрическая идентификация — определение параметров 
модели объекта. Проводится на базе математических моделей. 

Модель объекта представляет собой черный ящик, которая 
описывается зависимостью (функцией) выходного сигнала у или 77 от 
входного сигнала (рис. 1.1). = — случайные воздействия на объект, 
которые могут учитываться в модели. 







Процесс 





Рисунок 1.1 — Статистическая модель 


Количество входных сигналов не ограничено и зависит от 


математической сложности функции и учитывающихся в модели 
факторов. Входы могут называться факторами, входными сигналами, 
воздействующими сигналами, сигналами воздействия. 

Независимые входные сигналы х; в эксперименте называют 
факторами. В модель необходимо включить все существующие 
факторы, которые оказывают влияние на выходную величину у. 

Для отсеивания несущественных факторов используют метод 
случайного баланса. 

Факторы разделяются на количественные и качественные. К 
качественным факторам относятся разные вещества, разные 
технологические способы, аппараты, исполнители и т.д. 

При проведении эксперимента к факторам предъявляются 
следующие требования: 

• факторы должны быть управляемы, т.е. они могут быть 
установлены на определенный уровень и поддерживать уровень 
постоянным в течение всего опыта, в случае активного эксперимента. 
Если уровень фактора изменяется («плывет»), то следует использовать 
специальный метод, например, замер факторов по возможности 
должен быть более точным; степень точности устанавливается в 
зависимости от диапазона изменения факторов; 

• факторы должны быть однозначными. 

Количество выходов черного ящика — один, и может носить 
названия: выходной сигнал, функция отклика, реакция модели и, 
просто, выход. 

К выходному сигналу предъявляются следующие требования. 

1. Эта величина должна количественно оцениваться, т.е., 
например, поддаваться измерению. Если нет возможности 
количественно измерить его, то используют ранговый подход так 
называемое ранжирование. Величине п присваивают оценки, 
называемые рангами по заранее выбранной схеме, например 
двухбалльная, десятибалльная и т.п. В простейшем случае область 
определения „может содержать два значения: да и нет, хорошо и 
плохо (годная продукция – брак). 

2. Должна выполняться однозначность в статическом смысле. 
Это означает, что заданному набору значений х; должно 
соответствовать одно с точностью до ошибки экспериментальное 
значение 7. Если ошибка значительная, то необходимо повторять 
опыты. 

3. Должна оценивать эффективность функционирования 
системы. 

4. Желательно, чтобы выходная величина обладала свойством 
универсальности и полноты, т.е. всесторонне характеризовала объект; 


чтобы параметр у имел физический смысл, был простым и 
поддавался вычислению. 


1.4 Стратегия эксперимента 

Экспериментально-статистические методы в основном 
базируются на использовании пассивного и активного эксперимента. 

При пассивном эксперименте исследователь находится в роли 
пассивного наблюдателя. Эксперимент ведет сама природа. 
Экспериментатору приходится только фиксировать значения входных 
и выходных величин. Модели, полученные методом пассивного 
эксперимента почти не удается проверить на адекватность. 

При активном эксперименте исследователь вмешивается в 
процесс эксперимента путем варьирования уровней входных величин. 

В рамках активного эксперимента построение модели проходит 
следующие этапы: 

1) выбирается форма модели процесса; 

2) строится план эксперимента; 

3) проводится экспериментирование; 

4) дается анализ результатов эксперимента. 

На практике экспериментатору приходится чаще планировать не 
один, а несколько экспериментов, выполняя и анализируя каждый и, в 
соответствии с результатами, изменять план эксперимента. Стратегия 
такого эксперимента показана на рисунке 1.2. 






Экспериментирование 







План Иселрдувмый Анализ 
(реальный) 


эксперимента результатов 
процесс 





Форма 
модели 
процесса 







Рисунок 1.2 — Стратегия эксперимента 


В результате анализа результатов эксперимента может 
возникнуть необходимость исправления формы модели и плана 
эксперимента, тогда эксперименты повторяются вновь по указанной 
схеме. 
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1.4 Тактика эксперимента 

Стратегическое планирование применяется при любом виде 
экспериментов, как при натурных, так и машинных экспериментах. 
При машинных экспериментах дополнительно возникает 
необходимость обеспечения точности и достоверности результатов, а 
также количества прогонов. 

В связи с этим при тактическом планировании решаются 
следующие вопросы: 

1) определяются начальные условия и их влияние на 
достижение результатов процессов моделирования; 

2) обеспечения точности и достоверности результатов; 

3) уменьшение дисперсии оценок характеристики процесса 
функционирования систем; 

4) выбора правил автоматической остановки имитационного 
эксперимента. 

Рассмотрим каждый вопрос подробнее. 

Определение начальных условий усложняется тем, что 
имитационная модель работает эпизодически, т.е. только когда 
экспериментатор проводит эксперименты. Поэтому, начальный 
период работы искажается из-за влияния начальных условий запуска 
машины. Для решения этой проблемы исключают из-за рассмотрения 
информацию о модели, полученной в предыдущей части 
моделирования, либо начальные условия выбирают так, чтобы 
сократить время достижения установившегося режима. Эти приемы 
позволяют только уменьшить, но не свести к нулю время переходного 
процесса при проведении машинного эксперимента. 

Вопросы обеспечения точности и достоверности результатов, 
уменьшение дисперсии оценок характеристики процесса 
функционирования систем и выбора правил автоматической 
остановки имитационного эксперимента будут рассмотрены в 
разделах 5.3, 5.4, 4.5 соответственно. 


Контрольные вопросы 
Сформулируйте определение понятия «Моделирование». 
Сформулируйте определение понятия «Модель». 
Приведите общую классификацию моделей систем. 
Идентификация — это? 
Что такое фактор? 
Какие требования предъявляются к выходному сигналу? 
Сформулируйте понятие пассивного эксперимента. 
Сформулируйте понятие активный эксперимент 
Перечислите этапы построения модели. 
10. Какие вопросы решаются при тактическом планировании. 


они ыв- 
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2. Форма модели процесса 


2.1 Имитационное моделирование. Моделирующий алгоритм 

Для моделирования процесса на ЭВМ необходимо 
преобразовать математическую модель в специальный моделирующий 
алгоритм [8, 12, 14]. 

Существуют различные способы представления моделирующего 
алгоритма. 

Исторически первым из таких способов является запись 
алгоритмов при помощи операторных схем. Затем появились языки 
программирования, пакеты прикладных программ и специальные 
языки моделирования. 

Среди важнейших операторов моделирующих алгоритмов 
выделим следующие. 

1. Вычислительные операторы. Это В основном 
арифметические операторы. В операторных схемах моделирующих 
алгоритмов вычислительный оператор может описывать любую сколь 
угодно сложную и громоздкую группу операций. 

2. Операторы формирования реализаций случайных чисел, 
величин, событий и процессов (см. раздел 4.2, 4.3). 

Для имитации действия различных случайных факторов, 
сопровождающих исследуемый процесс, при моделировании 
возникает необходимость формировать реализации случайных 
событий, случайных величин и случайных функций. 

Для реализации всех этих перечисленных случайных групп 
необходимо иметь случайные числа. Операторы формирования 
реализаций случайных процессов решают задачу преобразования 
случайных чисел стандартного вида в реализации случайных 
процессов с заданными свойствами. 

3. Операторы формирования неслучайных величин. Они 
нужны при моделировании детерминированных математических 
моделей для реализации различных констант и неслучайных функций 
времени. 

2.1.1 Принципы построения моделирующих алгоритмов 


В большинстве имитационных моделей имитируется поведение 
системы на некотором отрезке времени. Поэтому при создании 
модели и выборе языка программирования важной задачей является 
определение механизма системного времени. Это необходимо для 
корректирования временной координаты состояния системы и для 
обеспечения согласованности различных блоков и событий в системе. 

Модель функционирует в искусственном времени и необходимо 
обеспечить появление событий в определенном порядке с 
надлежащими временными интервалами между ними. 
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2.1.2 Принцип фиксированного шага. Принцип переменного шага 


Существует два основных метода задания времени: с помощью 
фиксированных и переменных интервалов времени. Их иногда еще 
называют метод фиксированного шага 2 и шага до следующего 
события. По методу фиксированного временного шага отсчет 
системного времени ведется через заранее определенные временные 
интервалы. При использовании метода переменного шага состояние 
моделируемой системы обновляется с появлением следующего 
события. 

В непрерывных моделях используются механизмы 
фиксированных приращений временных интервалов. В большинстве 
моделей с дискретным изменением событий используется метод 
отсчета времени до следующего события. 

2.1.3 Принцип особых состояний 


Бывают случаи, когда соотношения математических моделей 
систем удается преобразовать таким образом, что можно определить 
последующее особое состояние системы по предыдущему особому 
состоянию или нескольким предыдущим состояниям. Для таких 
случаев моделирующий алгоритм может быть построен по так 
называемому принципу "особых состояний". Он отличается от 
принципа 2 только тем, что включает в себя процедуру определения 
момента времени, соответствующего следующему особому состоянию 
по известным характеристикам данного или предыдущего состояния. 

2.1.4 Принцип последовательной проводки заявок 


В системах массового обслуживания при моделировании 
процессов обработки заявок иногда удобно строить моделирующие 
алгоритмы по принципу воспроизведения истории отдельных заявок в 
порядке поступления их в систему. Алгоритм обращается к сведениям 
о других заявках лишь в том случае, если это необходимо для 
решения вопроса о дальнейшем порядке обслуживания данной заявки. 
Такие алгоритмы весьма экономны, однако они требуют весьма 
сложную логическую структуру. При моделировании систем 
массового обслуживания и дискретных производственных процессов, 
как правило, используют этот принцип, называемый иногда "принцип 
последовательной проводки заявок". Данный принцип, иногда в 
сочетании с принципом / , лежит в основе почти всех языков 
моделирования. 


2.2 Языки моделирования 
Для имитационного моделирования используются как 
универсальные, так и специальные языки [2,4]. Универсальные языки 
дают большие возможности программисту в смысле гибкости 
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разработки, отладки и использования модели. Однако это требует 
больших усилий, затрачиваемых на программирование. В этом случае 
усложняется порядок выполнения операций, отсчет системного 
времени и контроль. 

Специализированный язык отличается от универсального 
своими специфическими свойствами. 

К этим свойствам или требованиям можно отнести: 

• способность генерировать случайные числа; 

• возможность генерировать случайные величины и процессы; 

• возможность "продвигать" время либо на одну единицу 2, 
либо до следующего события; 

• способность накапливать выходные данные; 

• способность проводить статистический анализ 
накапливаемых данных; 

• способность распределять выходные данные по заранее 
заданным форматам; 

• способность выполнять идентификацию конкретных 
событий. 

Некоторые из языков имитационного моделирования являются и 
описательными языками. Они близки к естественному языку, поэтому 
имитационные модели, написанные на таком языке, легче 
воспринимаются руководителями и специалистами, не имеющими 
непосредственного отношения к программированию. 

Специальные языки цифрового имитационного моделирования 
делятся на две группы, соответствующие двум видам имитации: 
дискретных и непрерывных процессов. 

Для моделирования непрерывных процессов используются 
ЯЗЫКИ: 

• Динамо — для аппроксимации непрерывных процессов 
используют дифференциальные уравнения первого порядка; 

• (СМР, Міааѕ и др., которые помимо блочного построения, 
применяемого в языке Динамо, еще обладают мощью и удобством 
алгебраической и логической алгоритмизацией. 

Языки для моделирования дискретных процессов можно 
разбить на четыре категории: 

1) языки, ориентированные на действия (С51, Форсим ПУ и 
др.); 

2) языки, ориентированные на события (Симскрипт, Симком и 
др.); 

3) языки, ориентированные на процессы (Симула, 501); 

4) языки, ориентированные на потоки сообщений (СР55, 
ВО55). 
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В первой группе действия представлены в модельном времени 
как мгновенные. В этих языках нет регламентации действиям. Вместо 
этого исполнительные программы просматривают набор всех 
условий, от которых зависит появление какого-либо события. И 
только тогда, когда выполняются все контролируемые условия, 
происходит изменение состояния и сдвиг времени в этой части 
программы. 

Некоторые задачи удобнее программировать на языке 
ориентированном на события. При этом событие регламентировано. 
Регламентация обеспечивает наступление события именно в тот 
момент времени, когда динамическое состояние показывает, что 
сложились условия для его появления. 

Языки, ориентированные на процессы, объединяют достоинства 
первых двух языков, т.е. краткость языков, ориентированных на 
действия и эффективность языков, ориентированных на события. 
Написанная на этом языке программа работает как несколько 
независимых программ: одна посредством просмотра действий, 
другая посредством регламентирования событий. 

Что касается четвертой группы языков, то они на самом деле 
являются языками процессов и отличаются только схемным 
построением. 

Язык СР585 обычно выделяют в отдельную группу. Язык СР55 
представляет собой интерпретирующую языковую систему, которую 
применяют для описания пространственного движения объектов. 


2.3 Виды моделей 

Стратегия эффективного планирования эксперимента 
начинается с выбора формы модели процесса [7]. В зависимости от 
априорной информации о процессе различают две задачи: 

• построение модели при неизвестной структуре процесса; 

® нахождение параметров модели при заданной структуре. 

Априорная информация о процессе для первого случая почти 
отсутствует. При этом необходимо решать задачу выбора структуры 
модели процесса, устанавливать класс моделей, оценивать степень 
стационарности и линейности процесса и др. Решение подобных задач 
связано с большими трудностями. На заданном этапе приходится 
решать вопросы содержательного неформального эвристического 
характера. 

При известной функциональной зависимости между выходными 
и входными величинами задача является более простой и доступной. 
При этом форма модели известна. Экспериментатору остается 
построить план, провести эксперимент и сделать анализ результатов 
экспериментов. 
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2.3.1 Модели в статике 


Модели в статике бывают следующего вида (рис. 2.1). 





Рисунок 2.1 — Схема модели в статике 


1. Линейная по входным величинам х и по параметрам: 
п = Во + Вуху + В2х2 +...+ Вихи (2.1) 
2. Нелинейная по х и линейная по коэффициенту д: 
2 п 

п = Во + Вуху + В2х5 +...+ Вихя 

3. Линейная по х и нелинейная по коэффициенту р: 
2 п 

п = Во + Вуху + Вх. +...+ Вихи 

4. Нелинейная по х и пор: 


п= ей +6252 +. + епи 





, 
п= {Во + Врх + Водо +...+ ух, 

где 7] — реакция модели на входные величины Х. 

Встречаются модели и другого вида, например, 
тригонометрические, логические и др. 

Определение функциональной связи по опытным данным 
представляет трудную задачу и универсальных способов не 
существует. 

Можно указать простой способ определения вида моделей с 
одной независимой переменной: строят график полученных 
экспериментальных данных. Если экспериментальные точки образуют 
некоторую гладкую кривую, то форма этой кривой может подсказать 
аналитический вид кривой. 

Для проверки согласованности выбранной формы модели с 
экспериментальными данными преобразуют переменные к линейной 
форме (модели): 

у= &+ Вх 
Пример. Преобразовать переменные к линейной форме. 


1 1 
1. Модель 7 = 2 + 8 – . Принимаем х = —, у = 7. Тогда у= 0 + 2. 
"4 "а 


2 п=ау Вю п=юЮза+ Вову 
у= а + Вх;а =юга, х = 105 у 
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и Т.Д. 

Таким образом, параметры нелинейной модели можно 
определить по экспериментально статистическому исследованию 
модифицированной линейной модели. Однако, данный способ трудно 
применить в тех случаях, когда к выходной переменной модели 
добавляется ненаблюдаемая ошибка Є . 


Например, выражения п=ауй += и 1057 = 050+ Вю у+Е= 
неравнозначны, т.к. надо еще прологарифмировать ошибку & . 

При исследовании реального процесса всегда имеется 
возможность представить результаты экспериментальных данных 
множеством моделей. Чтобы выбрать одну из них, надо знать для 
какой цели находится модель, и какие требования к ней 
предъявляются. Очевидно, главным требованием следует считать 
способность модели предсказывать направление дальнейших опытов. 
Другим требованием является минимизация суммы квадратов 
отклонений наблюдаемых значений от вычисленных. Далее, 
желательно, чтобы преобразованная величина подчинялась 
нормальному распределению и чтобы дисперсия преобразованной 
случайной величины не зависела от самой величины. 

Важным требованием является простота модели. Если 
несколько различных моделей отвечают необходимым требованиям, 
то следует принять ту из них, которая имеет более простую 
зависимость. 

Рассмотрим графическую интерпретацию уравнения модели. 


^ у, 








> 
Хх 





Рисунок 2.2 — Графическая интерпретация уравнения модели 


На рисунке 2.2 изображена логарифмическая функция. На 
отрезке [хииХтх]она с определенной точностью может быть описана 
двумя уравнениями: 

У1 =Ю#рх (2.2) 
у? = 0х (2.3) 
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Уравнение (2.3) является линейным алгебраическим и его 
функциональная зависимость проще, чем уравнение (2.2), которое 
является трансцендентным. 

Поэтому, при прочих равных условиях более 
предпочтительными являются степенные ряды. Таким образом, выбор 
модели, как правило, следует искать среди алгебраических 
полиномов. Причем на первом этапе планирования эксперимента 
выбирают линейную модель в виде алгебраического полинома первой 
степени (2.1). 

2.3.2 Динамические модели 


Динамические модели можно подразделить на два больших 
класса: традиционные модели и современные модели или ещё их 
называют модели в пространстве состояний (рис. 2.3). 

К традиционным моделям относят простейшие инерционные 
модели, модели на базе передаточных функций и комплексных 
коэффициентов передачи. 


(0) 


у(і 2( 





х(0) 


— 


Детерминированная 
модель 











Рисунок 2.3 — Схема динамической модели 


Многие технологические процессы обладают инерционностью. 
Говорят, что такие процессы имеют «память». Инерционные модели 
бывают дискретные и непрерывные, линейные и нелинейные. 

Если модель линейная дискретная, то её представляют в виде 
суммы свертки (весовая функция): 

т-1 
УпТ]= У; е[т2Ат|: хпАЕ- (т + 1)Ат| + е[пЛ] (2.4) 
т=0 

где У[иТ]- дискретные значения выходного сигнала; 8[тАт]- 
дискретная весовая функция; ХИЛА – (т+1)А7|- дискретные 
значения входного сигнала; е[иД{- случайная помеха. 

В этой модели по существу нет параметров, она 
непараметрическая. Роль величин, которые необходимо определить из 
экспериментальных данных, играют значения ординат импульсной 
характеристики, которые рассматривают как коэффициенты 
регрессивной модели, а роль факторов здесь играют значения одной и 
той же входной величины, но в разные моменты времени. 
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Если линейная модель непрерывная, то модель будет типа 
интеграла свертки: 


[0,0] 
У(В) = [<()ха – т)ат +е(1) (2.5) 

0 
Модель (2.5) также как и (2.4) является непараметрической. Она 
содержит неизвестную весовую функцию е( т). На практике широко 
используется возможность представления весовой функции для 
стационарной системы в форме Релея-Ритца путем разложения 

функции в ряд по системе известных ортогональных функций: 


К 
807) = аф; (т,а) 
1=1 

где ф;(т,а) — заданная система базисных функций (фильтров), 
зависящая от параметра а. 

Такой прием делает модель параметрической. Теперь она 
содержит ограниченное множество параметров а, подлежащих 
определению. 

Нелинейные инерционные модели могут быть представлены в 
виде сумм или рядов Вольтерра. 

В импульсном варианте модель можно представить: 


У[ий = 80(иГ)+ )81(п,0):х(п- 0 + >, 5200,1, )-хт-д-хт-л 


і=0 і=0)=0 
(2.6) 
Для непрерывного объекта: 
[0.0] 
У(В) = 20(7) + [167 х(7-тр)ату + 
Т1 =0 
(2.7) 
[0.0] [0.0] 
+ | [82( 11,12). ха-т)-х@-то Ато +...+е(7) 
Т1 =072 =0 


В выражениях (2.6) и (2.7) ег и в соответственно весовые 
функции, их называют ядрами Вольтерра первого и второго порядка; 
80 - составляющая не связанная с входным сигналом. 

Непосредственное определение ядер по опытным данным 
представляет собой сложную задачу. Поэтому их обычно 
аппроксимируют путем разложения в ряд по системе известных 
ортогональных функций: 

К 
8171,4) = У,а;ф;(ту,9) 
1=1 


19 


К К 
520,71,72,9) = а ФЕ, @)-ф (2,0) 
Е] 

Теперь задача построения модели сводится к определению 
параметров весовой функции по косвенным экспериментальным 
данным. 

Модели на базе передаточных функций 

Для непрерывных динамических объектов выходная величина 
определяется интегралом свертки (2.5). 

Применим преобразование Лапласа. На основании теоремы 
свертывания имеем: 


у(р) = 8 (р) х(р)+е(р) (2.8) 

Для импульсного динамического объекта, согласно выражения 
(4.4) после 7- преобразования получаем: 

у(х) = 8(<) - хе) + е(х) (2.9) 

Модели (2.8) и (2.9) имеют вид регрессионных зависимостей. 
Параметры этих моделей, которые необходимо определить из 
экспериментальных данных, содержатся в выражениях передаточных 
функций. 

Модели на основе комплексного коэффициента передачи 

К модели частотного вида можно перейти путем замены в 
выражении (2.8) аргумента р на ју: 

УС Ли) = 8( 1) - хр) +е м) (2.10) 

В роли аргументов в моделях (2.8), (2.9), (2.10) выступает не 
время, а соответственно параметры преобразования: г, р, ји. Все эти 
модели линейны по входным сигналам, но, как правило, не линейны 
по параметрам. 

Современные модели (модели в пространстве состояний) 

В настоящее время все чаще переходят к математическому 
описанию динамических систем в описание систем в пространстве 
состояний. При этом используют конечно — разностные и 
дифференциальные уравнения в форме Коши, т.е. разрешенных 
относительно первых разностей и первых производных. Описание 
систем в пространстве состояний позволяет с единых позиций 
рассматривать различные системы: линейные, нелинейные, 
дискретные и непрерывные. 

Модели в виде конечно — разностных уравнений 

Пусть динамический объект описывается линейным конечно- 
разностным уравнением 1-го порядка с постоянными 
коэффициентами. Далее считаем, что наблюдения производятся в 
дискретные равноотстоящие моменты времени. Тогда модель можно 
представить выражением (2.4). 

Примером таких моделей являются дискретно — непрерывные 
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модели. 
Модели в виде обыкновенных дифференциальных уравнений 
Обыкновенное линейное дифференциальное уравнение /-го 
порядка с постоянными коэффициентами имеет вид: 
а! у 41у 42у аїх 49! 
со Г 6] 11 Са. === 
аї аї 








(2.11) 
Уравнением (2.11) описываются рассмотренные выше 
непрерывно- детерминированные модели (О- схемы). 


Контрольные вопросы 

1. Какие способы представления моделирующего алгоритма 
Вы знаете? 

2. В чем заключаются принципы фиксированного и 
переменного шага? 

3. В чем заключается принцип особых состояний? 

4. В чем заключается принцип последовательной проводки 
заявок? 

5. Перечислите отличия специализированных языков от 
универсальных языков. 

6. Какие виды моделей в статике Вы знаете? Перечислите их 
особенности. 

7. Какие виды динамических моделей Вы знаете? Перечислите 
их особенности. 
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3. План эксперимента 
Будем считать, что математическая модель выбрана и 
представляет собой зависимость 
п = Во + Вх + 5х2 +...8,х, 
Выходная величина ту может иметь и более сложную форму 
переменных х,. 
Если поставлена задача найти оптимальные условия для и, то 


выходную величину называют параметром оптимизации. 

В качестве выходного параметра 77 может быть температура, 
давление или более сложный показатель качества выпускаемой 
продукции. 

В тоже время в планировании могут участвовать сложные 
факторы, например: соотношения между компонентами; 
соотношения, определяемые алгоритмом функционирования и т.д. 

С помощью сложных факторов можно учесть некоторые 
динамические свойства объекта, исследуя его в статическом режиме. 
Например, если какое-то измерение во времени имеет сложный вид 
кривой, то целесообразно в качестве фактора обозначить номер 
кривой. Тогда различные варианты кривых будут рассматриваться как 
уровни. Таким образом можем представить сложный фактор — 
функцию с помощью простых однозначных факторов. 

Выбор экспериментальной области факторного пространства 
связан с тщательным анализом априорной информации. В этой 
области находят локальную подобласть для планирования 
эксперимента. 

При наличии априорных сведений о процессе за основной 
уровень принимаются те значения факторов, при которых выходная 
величина принимает лучшее значение. Построение плана 
эксперимента сводится к выбору экспериментальных точек, 
симметричных относительно основного уровня. 


3.1 Выбор интервалов варьирования 

После выбора основного уровня выбирают интервал 
варьирования для каждого из факторов. 

Интервалом варьирования факторов называют некоторое 
число, прибавление которого к основному уровню дает верхний, а 
вычитание — нижний уровни факторов. То есть это расстояние на 
координатной оси между основным и верхним (либо нижним) 
уровнем. 

Для упрощения записи условий эксперимента и обработки 
экспериментальных данных масштабы по осям выбираются таким 
образом, чтобы верхний уровень соответствовал +1, а нижний -1, 
основной — нулю. Уровни, записанные в таких обозначениях, 


22 


называются кодированными р 


с (3.1) 
х =“ 
Ј 
А, 
где х, — кодированное значение фактора; / — номер фактора; 
х. натуральное значение фактора; х. — натуральное значение 


Ј ЈО 
основного уровня; Я; —— интервал варьирования. 


Пример. Дано интервал варьирования 14;=2; натуральное 
значение основного уровня х;,=3 (рис. 3.1) 





Рисунок 3.1 — Кодирование уровней 


Е 5 5-3 
Верхний кодированный уровень: х, = > 1; нижний 
В 1-3 
кодированный уровень: х, = 25. =-1. 
Экспериментальные планы, в которых все факторы выбираются 
К 
на двух уровнях, называются планами типа 2°, где К — число 


факторов (рис. 3.1). Встречаются планы 3“ и более числа уровней. 

На выбор интервала варьирования накладываются естественные 
ограничения сверху и снизу. Интервал варьирования не может быть 
меньше ошибки экспериментатора, которую он допускает при 
фиксации уровня фактора. В противном случае верхний и нижний 
уровни будут неразличимы. С другой стороны интервал варьирования 
ограничен сверху пределами области определения. 

После определения интервалов варьирования выбирают 
значения факторов, т.е. число уровней для каждого фактора. Число 
уровней всегда ограничено снизу и не может быть меньше двух, т.е. 
если принять один уровень, то фактор окажется постоянным во всех 
опытах. 

Рандомизацией называется любая процедура, обеспечивающая 
случайный порядок проведения опытов. Такую процедуру 
осуществляют с помощью случайного выбора номера опыта. При 
проведении эксперимента всегда следует применять рандомизацию. 


3.2 Полный факторный эксперимент 
Эксперимент, в котором реализуются все возможные сочетания 
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уровней факторов, называют полным факторным экспериментом 
К 

(ПФЭ). При двух уровнях имеем ПФЭ типа 2`. Число опытов для 

данного случая будет равно 


К 
м2 
Условие эксперимента записываются в виде таблицы. Строки её 
соответствуют различным опытам (вектор — строка), столбцы — 


значениям факторов в кодированном виде (вектор — столбцы). Такие 
таблицы называются матрицами планирования эксперимента 
(МПЭ). 
Составим МПЭ для двумерной модели 
у = Вх + 5х5 +Е (3.2) 
на двух уровнях 2? (табл.3.1). Число опытов №=2?=4. 
Таблица 3.1 





План эксперимента можно представить геометрически (рис. 

2 У 

3.2). Для плана 2“ каждая комбинация факторов представляет собой 
вершину квадрата. 


2тп 





0 


. Х Хх 
1тіп 10 НЕХ 1 


Рисунок 3.2 — Геометрическое представление ПФЭ 


В области определения факторов находят точку 
соответствующую основному уровню. Через эту точку проводят 
новые оси координат, параллельно осям натуральных значений 
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факторов. Затем выбирают масштабы по новым осям для каждого 
фактора согласно выражению (3.1). 

В МПЭ вводится фиктивный столбец хо для учета свободного 
члена Бо. 


Коэффициенты Ду, В, 3, оцениваются согласно выражений: 
Ху, Ху, ху, 
Бо Во р. = 
М М М 
3.2.1 Свойства полного факторного эксперимента 2“ 








К свойствам МПЭ относятся те, которые определяют качество 
модели, т.е. эти свойства делают оценки коэффициентов модели 
наилучшими. Первые два свойства вытекают из построения матрицы. 

1. Симметричность относительно центра эксперимента. 


Алгебраическая сумма элементов столбца каждого фактора равно 
м 

нулю Ух; =0, где /- номер фактора, № - число опытов. 
=] 
2. Условие нормировки. Сумма квадратов элементов каждого 


м 
столбца равна числу опытов У х п = №. 
11 
3. Ортогональность матрицы. Сумма почленных произведений 


М 
любых двух векторов — столбцов матрицы равна нулю У х ихш = 0, 
= 
гле: 50у), 
Ортогональные планы делают эксперимент более эффективным. 
Ортогональность плана позволяет получить оценки для 
коэффициентов уравнения регрессии независимые друг от друга. 
Иными словами ортогональность характеризует отсутствие 
корреляции между факторами. Однако если имеет место 
нелинейность, то столбцы взаимодействий окажутся неразличимы, 
закоррелироваными с некоторыми столбцами линейных эффектов. 
Это приводит к тому, что по результатам данного эксперимента 
становится невозможным разделить коэффициенты регрессии между 
линейными и нелинейными факторами. 
4. Рототабельность планов. Это такие планы, для которых 
дисперсия У одинакова для всех точек пространства переменных х, 


лежащих на одинаковых расстояниях от центра. 
3.2.2 Выбор модели при проведении полного факторного 
эксперимента 
Планируя эксперимент на первом этапе, всегда стремятся 
получить линейную модель. Для двух факторов модель представляют 
в виде выражения (3.2). Однако не всегда экспериментатор имеет 
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гарантии, что в выбранных интервалах варьирования процесс 
описывается линейной моделью. Часто встречающийся вид 
нелинейности связан с эффектом взаимодействия между факторами. 
ПФЭ позволяет оценить кроме коэффициентов при линейных 
эффектах коэффициенты взаимодействия. Для этого перемножают 
соответствующие столбцы. Тогда уравнение принимает вид 
у = Вохо + Вх + Вьх> + В1охихо += (3.3) 
МПЭ с учетом фактора взаимодействия для ПФЭ 2? показана в 
таблице 3.2. 
Таблица 3.2 


МПЭ для ПФЭ 2? 





Коэффициенты уравнений регрессии (3.2, 3.3) оцениваются 
следующим образом: 
М 


2.5, 2 у рат 


а ‚В, >Ь, У ан аа. 





По столбцам х; и х› осуществляют планирование, что же 
касается столбцов , хои хух› ‚то они служат только для расчета. 

При составлении МПЭ руководствуются следующими 
правилами: 

• располагают, если имеется соответствующая информация, 
факторы в матрице в порядке убывания степени их влияния на целую 
функцию; 

® стремятся ВЫПОЛНИТЬ требования рандомизации 
варьирования уровней. 

При составлении матрицы уменьшают частоту чередования 
уровней при переходе от ху к 2, от х2кхз и т.д. каждый раз вдвое. 

Рассмотрим пример составления МПЭ для трех факторного 
полного эксперимента. В качестве уравнения регрессии берем 
неполную квадратичную модель. 

3 3 
р р Р (3.4) 
і=1 і<ј 


Введем обозначение переменных х через <, тогда 


< 
ӯ = УЬ; (3.5) 
і=1 
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3 3 3 
где рУ В: = БиХ. Бе, 
і=1 і=1 і<ј 
Составим МПЭ. № = 2° = 8 (табл. 3.3). 
Таблица 3.3 
МПЭ для ПФЭ 2° 
Номер Кикини хІх; ххх 
опыта 


0 





сав 
раги Е 
а Еа 


После экспериментирования проверяют однородность 
дисперсии, значимость коэффициентов и адекватность модели. 


3.3 Дробный факторный эксперимент 

Во многих реальных процессах некоторые факторы 
взаимодействия могут отсутствовать. И тогда ПФЭ будет обладать 
избыточностью опытов. 

Рассмотрим пути минимизации числа опытов. 

Обратимся к уравнению (3.3). Если располагаем сведениями о 
том, что в выбранных интервалах варьирования процесс в статике 
может быть описан линейной моделью, то достаточно определить три 
коэффициента Бо ‚Р; ,Р›. В результате остается одна степень свободы, 
т.к. имеем четыре опыта, а количество констант три. Используем эту 
степень свободы для минимизации числа опытов. При линейном 
приближении Р; 0 и тогда вектор — столбец хух› может быть 
использован для нового фактора х;. 

При этом эксперименте появляются смешанные оценки 

р >В + 853,0 Э В> + Вз,Ьз > Вз + Ву, (3.6) 
т.е. столбцы. 


Таблица 3.4 


МПЭ для ДФЭ 2“! 





Пример. Допустим х; и х›х; между собой неразличимы. Однако 
парные взаимодействия в линейной модели незначительны. Зато 
вместо восьми опытов для изучения влияния трех факторов можно 
поставить только четыре опыта, т.е. вместо ПФЭ 2° получаем 2°. В 
теории эксперимента 2° называют полу —репликой. В общем случае 
имеют дело с дробной репликой. А факторный эксперимент называют 
дробным (ДФЭ). 

Для уяснения принципа составления МПЭ ДФЭ введено 
понятие определяющего контраста. Он позволяет определить, какие 
оценки смешаны друг с другом, не изучая МПЭ для выявления 
совпадающих столбцов. Для этого берут символичное обозначение 
произведения столбцов равного +1 или -1. Это и называют 
контрастом. Чтобы определить какой эффект смешан с данным, 
нужно помножить обе части определяющего контраста на столбец, 
соответствующий данному эффекту. 

Пример. Пусть имеем три фактора х; ,х› ‚хз. При построении 
полуреплики 2° имеется только две возможности приравнять хз либо 
к «+хих2», либо к «-х1х2» (табл. 3.5). 

Таблица 3.5 


МПЭ для ДФЭ 2°! 





Возьмем в качестве определяющего контраста —1= хух›х; 


„Гогда – ху = а Учитывая, что х2 =1 получаем ху = –х›х3. 
Теперь возьмем за определяющий контраст +1 = хух›23. 
Получаем: ху = х›х3,х5 = хухз, хз = хх. Эти выражения показывают, 
что коэффициенты линейного уравнения будут оценками (3.6). 
Соотношение, показывающее с какими из эффектов смешан 
данный эффект, называется генерирующим соотношением. 


4-1 
При выборе полуреплики 2° возможны восемь генерирующих 
соотношений: 


1.20 ЕЕ ЕО 

2. ха = хх 6. х = хх; 
ЕО: ГЕ асвет 
4. ха = хо хз Е 


Разрешающая способность этих полуреплик различна. Реплики 
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1 — 6 имеют по три фактора и носят название планов с 
расширяющей способностью Ш ( по наибольшему числу факторов в 
определяющем контрасте). Реплики 7 — 8 имеют по четыре фактора и 
обладают максимальной разрешающей способностью. Их называют 
главными репликами. Всегда стремятся выбрать реплику с 
наибольшей разрешающей способностью, т.к. чем больше эффектов 
взаимосвязано, тем точнее окажется полученная модель. 

Однако, если имеется информация об эффектах взаимодействия, 
то реплики нужно выбирать с ее учетом. 

Реализация МПЭ ДФЭ ничем не отличается от реализации МПЭ 
ПФЭ. Методика оценки значимости коэффициентов и проверка 
адекватности модели проводится также как и в ПФЭ. 


3.4 Обобщающий определяющий контраст 
Рассмотрим на примере исследование модели с пятью 
факторами. Возьмём реплику 27”. Получаем 8 опытов вместо 32. 
Возможны 12 решений, если приравнять х; парному 
взаимодействию, а х; — тройному: 


1. Ха=Х1Хэә Х5=Х1Х2Х3 

2, Ха=Х1Хэә Х5=-Х1Х2Х3 

3. Ха=-Х1Х2 Х5=Х1Х2Х3 

4. Ха=-Х1Х2 Х5=-Х1Х2Х3 

5. Х4—ХХз Х5=Х1Х2Х3 

б. Х4—Х1Х3 Х5—-Х1Х2Хз 

7. Ха=-Х1Х3 Х5=Х1Х2Х3 

8. Ха=-Х1Х3 Х5=-Х1ХэХ3 

9. Ха—Х2Х3 Х5=Х1Х2Х3 

10. Ха—Х2Х3 Х5=-Х1ХэХ3 

11. Х4=-Х2Хз Х5=Х1Х2Х3 

12. Х4=-Х2Хз Х5—-Х1Х2Хз 
Допустим выбран первый вариант. Тогда определяющими 
контрастами будут: 1=хлхх. , 1=хдохх5. Перемножим эти 


определяющие константы, получим третье соотношение: /=хзхХ5. Для 
того чтобы полностью охарактеризовать разрешающую способность 
реплики, вводят понятие обобщающего определяющего контраста: 
1= хах2х4= ХЗХАХ5ЕХИХ2Х 5. 
Система смешивания столбца определяется умножением 
обобщающего определяющего контраста последовательно на хр,х»,Хз: 
мА ХХХ ХХХ 5, 
ХХ МХА МХ3Х 5, 
Хз = ХІХХЗХА Ла 1Х2Х 5, 
ХХХ Хзх5 М2 34 5, 
Х5-ХАХХАХ5 Хз ХХХ 3, 
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ХІХ Ха ХХәХЗХАХ5 ХХ... 

Если при выбранной реплике некоторые коэффициенты 

получаются отличными от нуля, например: 

Бо —> Ву + Ва + 235 + Воздѕ 
то ставят вторую серию опытов с другой репликой, например, берут 
вариант 4. 

Дробные реплики находят широкое применение при получении 
линейных моделей, причем, целесообразность применения их 
возрастает с ростом количества факторов. Эффективность применения 
дробных реплик зависит от выбора системы смешивания линейных 
эффектов с эффектами взаимодействия. 


3.5 Планирование экспериментов при построении полной 
квадратичной модели 

В уравнениях (3.2), (3.3), (3.4) учитывались только линейные 
эффекты и эффекты взаимодействия. В некоторых случаях 
существенными могут оказаться коэффициенты при квадратных 
переменных, их кубов и т.д. 

Для двухфакторного эксперимента модель может быть 
представлена выражением 


у = рохо + руху + Рх) Е Вох, х> + Вх + Ь.>х2 (3.7) 
Полученные вектор - столбцы хѓи х2 являются единичными 
столбцами, совпадающие друг с другом и с фиктивным столбцом ху. 


Эти столбцы неразличимы, поэтому нельзя сказать за счет чего 
получилась величина Р,. Очевидно, она включает в себя значения 


свободного члена Во и вклады квадратичных членов. Символически 
это можно записать: 


к 
в э В+ У, В,; 


і=1 
Для квадратичной модели получается следующая система 
смешивания: 
ро ә Во + Ви, +2, ә Ву, 5,65 Вә. 
Следовательно, планирование эксперимента на двух уровнях не 
дает возможности получить раздельные оценки коэффициентов при 
квадратичных членах и фиктивной переменной х,. 


Согласно теории интерполяции, для решения задачи 
нахождения раздельных оценок число уровней каждой из 
независимых переменных должно быть на единицу больше степени 
интерполяционного полинома. Для полинома второй степени число 
уровней должно быть равно трем. 

Однако применение методов ПФЭ плана 3" не является 
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рациональным из-за резкого увеличения опытов эксперимента. 
Поэтому разработаны специальные методы построения планов 
второго порядка. 

Например, в качестве двухфакторных планов второго порядка 
могут служить планы, представляемые вершинами и, по крайней мере, 
одной центральной точкой любого (п-[) мерного правильного 
многоугольника (который можно вписать в круг). 

Пример. Имеем восьмиугольный план (рис. 3.3, табл.3.6). 

Этот пример можно обобщить на случай получения планов 
второго порядка. Для этого к ПФЭ типа 2" добавляется центральная 
точка с координатами (0,0....0) и, так называемые, звёздные точки с 
координатами (0,0,..., + с,...,0), лежащие на сфере диаметра 2а. Т.е. 
план ПФЭ достраивается до плана второго порядка. Такой план 
называется композиционным планом. 





Рисунок 3.3 — Восьмиугольный план эксперимента 


Таблица 3.6 
Матрица композиционного плана двухфакторной модели 


7 
План 2 
представлен 


квадратом 
АВСД 


План 
представлен 
звёздными точками 
ММКГ, 
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Для трехфакторной модели МПЭ будет иметь вид, 
представленный в таблице 3.7. 


Таблица 3.7 
Матрица композиционного плана двухфакторной модели 


+|+|+|+|+|+ |+ 
к к к | 


реч чучу чуча ГЕ 
+ 
— 


Ганя 
НЕК 
ИЯ 
ЕЮ 
1-28 
ЕВ 
ВЕЕ 
ЕЕ 
о 
ЕВ 
го) 
28 





Добавление двух сфер, образованных звездными точками и 
центральной точкой, к ПФЭ позволяет получить раздельные оценки Бо 
и Б. Все три сферы образуют композиционный план второго порядка. 

В зависимости от критерия оптимальности плана, различают 
ортогональное композиционное планирование и рототабельное 
композиционное планирование. 

План, приведенный в таблице 3.6, является рототабельным и 
обеспечивает получение раздельных оценок Рои В, 

3.5.1 Ортогональное центральное композиционное планирование 


Критерием оптимальности является ортогональность столбцов 
матрицы планирования. В силу этого свойства все коэффициенты 
модели определяются независимо друг от друга. 

Ортогональность столбцов хо и ох’ достигается путем 
преобразования вида: 


д =1 0 
32 =52- 1 =х2 52 (3.8) 





С учетом выражения (3.8) условие ортогональности 
выполняется: 








И М 

Из условия ортогональности выбирают координату © звездной 
точки для п — независимых переменных. 

Таблица 3.8 соответствует ортогональному композиционному 
плану для количества переменных и=2,3,4 , где №, — число звездных 
точек; № — число точек в центре эксперимента; № — количество 
точек куба (гиперкуба) при ПФЭ; № — общее число точек факторного 
пространства. 


Таблица 3.8 
Количество точек ортогонального композиционного плана 





Составим матрицу ортогонального планирования для 
трехфакторного эксперимента. 
В таблице 3.9 для точек ПФЭ 


Таблица 3.9 
Матрица ортогонального планирования для трехфакторного 
эксперимента 
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Где 1-8 опыт соответствует ПФЭ 2°; 9-14 опыт соответствует 
звездным точкам; 15 опыт - центральная точка. 

Уравнение регрессии для полного квадратичного полинома 
будет иметь вид: 


у= +УБх, + Ухх, +) = 
1=1 і=1 


У 


=, + хор + УВ хх, + ХА. (3.9) 


1,1—1 
* и —2 
где 2, = 0, – Ва г, п- количество факторов. 
і=1 
ж 
Дисперсия коэффициента Бу будет оцениваться по формуле: 
2 _ а2 2 —2 

5 = 5, + 2.5, ХЕ (3.10) 

Коэффициенты Ё;, Би Б; определяются по формулам: 
Е М = 2 № = 
ХНУ, Хх иі Уи Хау 

Ь, =2—__— ла а | р; Ееее 


ху 


=1 


9 


=. 


Р 


р 


— ра 


где и=1, № - число опытов; у, = ‚ Р- число прогонов. 


Р 

Из-за неодинаковой дисперсии коэффициентов регрессии 
критерий ортогональности является недостаточно сильным критерием 
оптимальности для планирования второго порядка. Поэтому точность 
предсказания выходной величины в различных направлениях 
факторного пространства неодинакова. 

Лучшим методом планирования является такой метод, который 
обеспечивает одинаковую точность во всех направлениях на 
одинаковом расстоянии от центра. Таким методом является 
рототабельное композиционное планирование. 
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3.5.2 Рототабельное композиционное планирование 


Критерием оптимальности в рототабельном планировании 
является условие с; =сопѕїі при одинаковом удалении точек 


эксперимента от центра, т.е. А = сопзі. 

Если имеются двухфакторные планы, то, как уже было 
отмечено, типичными примерами рототабельных планов являются 
планы, представляемые вершинами и, по крайней мере, одной 
центральной точкой любого (п-1) — мерного правильного 
многоугольника, который можно вписать в круг (рис. 3.3). 

Композиционные центральные рототабельные планы также как 
и ортогональные состоят из трех сфер: сфера нулевого радиуса - 
центральные точки; сфера точек куба или гиперкуба и сфера звездных 
точек. Равномерность расположения точек на сфере приводит к 
вырожденным матрицам. Для устранения вырожденности используют 
сферу нулевого радиуса с несколькими центральными точками. 


Таблица 3.10 
Количество точек рототабельного композиционного плана 





где № - число звездных точек; № - число точек в центре 
эксперимента; №. - количество точек куба (гиперкуба); № - общее 
число точек факторного пространства. 

Матрица планирования рототабельного плана второго порядка 
для трехфакторного эксперимента будет представлена в таблице 3.11. 


Таблица 3.11 
Матрица рототабельного планирования для трехфакторного 
эксперимента 
































Эксперимент проводится аналогично ПФЭ, однако оценки 
коэффициентов а по своим ов 


м 0002) - -2асў хия Р (3.11) 
і=1 ј=1 


СА = 
Ь,; =. С2[(0 + 2)4 – уузу кс! (1 дўи х. -21СУ хоў; 
і=1 ]=1 Ј=1 











(3.12) 
> 
ра Уу; (3.13) 
М1 Ле Л 
с 
Ь.. = У, (3.14) 
4 М = Д 
К 
У А пм > М№,Р,, 
С=- СУТ ТЛ а "= (3.15+3.17) 
п+ —п 
я (п+2)(У МР2)? 
Ј=і1 у=] 


где М, - число точек на сфере радиусарР,,; к - число сфер (к=3). 

Рототабельные планы нашли широкое применение на практике. 
Однако с точки зрения математиков, занимающихся развитием 
математической статистики, выбор такого критерия представляется 
мало обоснованным. Он не вытекает логически из тех идей, на 
которые базируется математика. Например, как выбрать расстояние о 
до звездных точек; не все пространства независимых переменных, 
отведенное для эксперимента, используется в композиционных 
планах и т.д. Невозможно из множества рототабельных планов при 
одном и том же числе факторов выбрать лучший, т.к. не было 
критерия оценки. Поэтому этот критерий стали относить к эмпирико- 
интуитивным критериям. 
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Наряду с развитием планирования экспериментов, основанных 
на эмпирико-интуитивным критериях Бокса в США стало 
развиваться второе чисто теоретическое направление, которое 
связывают с планом ученого Кифера. Он установил связь между 
некоторыми критериями оптимальности; теоретически доказал, что 
для отдельных видов регрессии одни и те же планы могут отвечать 
сразу нескольким критериям оптимальности. 


3.6 Разбиение матрицы планирования 2“ на блоки 

При проведении серии экспериментов, осуществляемых в 
течение некоторого периода времени или с использованием 
нескольких партий веществ может произойти изменения условий 
эксперимента. Для предотвращения смещения среднего уровня 
отклика точки экспериментального плана рекомендуется выбирать в 
случайной последовательности. Однако в некоторых случаях влияние 
внешних условий можно избежать путем подходящего разбиения 
матрицы на блоки по сырью или последовательности испытаний. Так, 
при наличии двух партий сырья матрицу 2° можно разбить на два 
блока таким образом, чтобы эффект сырья сказался на величине 
трехфакторного взаимодействия х/хохз (табл. 3.12). В этом случае все 


линейные коэффициенты и парные взаимодействия будут 
освобождены от влияния неоднородности сырья. 


Таблица 3.12 


МПЭ ПФЭ 2° 


арры әрәр 
ЕСЕ Е ЕИ ИЕ ЕЕ 
ЕЕЕ ЕА Е ЕЕЕ ВЕЕ ЕЗ 
Е ааа ЕЕ] 





Уз += 





Различие в сырье рассматривают как новый фактор х4. 





Матрица 2°, представленная двумя блоками: первый блок 
(ххохз = +1) и второй блок (хухохз = -1), является полурепликой 
2*' с определяющим контрастом 1 = х|хохзхд. 

Коэффициенты модели будут равны 
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1 
Бо и +=)+(у2+=)+(у3+=)+(у4+=)+ у5 + уб + у7 + 53| 
с 
В Во "5 


Ы НС +=)+(у2+=)-(у3+=)+(у4+=)- У5 + уб – У7+ 53| 


Ь> —> 2 
Ьз > 3 
Ы2 —> 12 
руз —> Віз 
Боз —> Воз 


В 23 = +=)+(у2+=)+(уз3+=)+(у4+=)- У5 – уб – Ут – 8] 


& 
Буэз >> 0123 2 
Эффект сырья отразился на подсчете свободного члена ру и 


коэффициента рэз. 

Аналогично можно разбить на два блока любой эксперимент 2". 
Нужно только выбрать такой элемент взаимодействия, которым 
можно безболезненно пожертвовать. Обычно, при отсутствии 
априорных сведений выбирают взаимодействие самого высокого 
порядка х] Х2х3 для О. ‚ ххохзхд ДЛЯ 2 хүхэхзх4х5 ДЛЯ х5 ИТ.Д. 

Если имеется четыре источника неоднородности, которые могут 
заметно исказить результаты эксперимента, то матрицу 2“ разбивают 
на четыре блока таким образом, чтобы линейные эффекты были 
освобождены от влияния межблокового эффекта. 

Схемы разбиения на блоки рототабельных планов второго 
порядка так же нужно выполнить с учетом того, чтобы различие 
между блоками не влияло на коэффициенты при линейных, 
перекрестных и квадратичных членах. Для того, чтобы не было 
искажения результатов эксперимента и результаты не зависели от 
блоковых эффектов число центральных точек должно быть выбрано 
подходящим образом. 


3.7 Критерии оптимальности планов 

Минимизацию чувствительности математической модели 
процесса к случайным факторам проведем с учетом дисперсионных 
оценок. 

Предположим, что вид модели известен. Требуется найти 
лучшие оценки коэффициентов уравнения регрессии. 

Если это уравнение линейное относительно параметров, то 
оценки уравнений регрессий определяются из системы нормативных 
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уравнений 
ХТУИХЬ = ХГУ, 
= (УГ 
где Х = {х н) — матрица независимых переменных; Х` = {х я) 
— матрица, получаемая транспонированием матрицы Х ; Ү = ХБ — 


уравнение регрессии в матричной форме; ИУ — функция веса; 
и тогда оценки коэффициентов 
р= (хх ) (ху). 

Будем предполагать, что дисперсия выходной величины У не 
зависит от независимых переменных Х и поэтому примем функцию 
веса И = 1. 

С учетом общего вида функциональной зависимости Ү от Х 
получаем следующие оценки: 


Ь= (ЕТ) Ету ий (3.18) 
ЛЕ Л2 .. Л Л 
Ри: 120. 5. Ри |0 


где Ё = = 5), убх) — известные функции 


Хп 2 52. Тик т 


входных переменных ху, х›,...,х,, Т.е. 











к 
т = УЬ; (х), (3.19) 
а 


= (бо... )- Обозначим 7 = ЕТү, 


Если Е - транспонированная матрица Е, то 
1 
с- (ЕТ) (3.20) 


—— ковариационная матрица. Для линейной модели С = (х УХ с 
В результате действия случайных факторов экспериментатор 
получает математическую модель отличную от (3.19): У, =; + є;. 
Определим чувствительность решения (3.18) к небольшим 
изменениям выходной величины У, а следовательно и 7 за счет 
случайного фактора: 


эр = = С (3.21) 
Из выражения (3.21) следует, что чувствительность решения "Р" 
по отношению к случайным факторам выражается ковариационной 


матрицей С. 
План эксперимента выбирается с учетом минимизации 
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чувствительности решения (3.18) по отношению к случайным 
воздействиям. Однако в том виде, в каком она предоставлена в 
выражении (3.21) не представляется возможным дать количественную 
оценку чувствительности для разных планов Е. Поэтому в качестве 
числовой характеристики чувствительности принимают некоторую из 
числовых характеристик ковариационной матрицы С. В зависимости 
от этого получают ряд критериев оптимальности. 

1. План называется А — оптимальным, если он минимизирует 
сумму квадратов главных полуосей эллипсоида рассеяния оценок. 
Известно, что при нормальном законе распределения результатов 
наблюдений, поверхностями постоянных значений функции 


р 


ПЛОТНОСТИ распределения Р| — ЯВЛЯЮТСЯ ЭЛЛИПСОИДЫ. 


Характеристики эллипсоида полностью определяются элементами 
матрицы С. Для двухфакторного эксперимента дадим 
геометрическую интерпретацию (рис. 3.4). 








В Ь, 


Рисунок 3.4 — Плоскость постоянных значений функции плотности 
при нормальном законе распределения 


Через 5 ,С обозначена сумма диагональных элементов матрицы 


К 
5С = У АС. 
і=1 
На языке матричной алгебры минимизация суммы квадратов 
главных полуосей эллипсоида означает минимизацию следа (суммы 
диагональных элементов) ковариационной матрицы уравнения 
регрессии. Еще можно встретить название минимизация средней 
дисперсии. 
2. План называется О — оптимальным, если он минимизирует 
величину определителя матрицы С. Это соответствует минимизации 
обобщенной дисперсии. 
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3. Е — оптимальный план минимизирует максимальную ось 
эллипсоида рассеяния. Или еще можно назвать минимизацией 
максимального собственного числа матрицы С или минимизации 
максимального характеристического значения ковариационной 
матрицы уравнения регрессии. 

4. © — оптимальный план обеспечивает наименьшую по всем 
планам максимальную величину дисперсии предсказанных значений 
функции отклика в области факторного эксперимента. Эти планы 
называют минимаксными. 

С — план можно определить из выражения 


а) "(ДЕТ Е)" (к)? (3.20) 
План будет С-оптимальным, если 
а) тип тах (ОДЕР) ђе? (3.23) 
Определитель матрицы С является наиболее полной числовой 
характеристикой чувствительности. Г — оптимальным план 
учитывает все элементы матрицы С. Поэтому он лучше Аи Е — 
оптимальных планов. Что касается С — оптимальным, то он 
эквивалентен О — оптимальному плану. р — оптимальный план 
минимизирует объем эллипсоида рассеяния оценок коэффициентов 
модели. 
Таким образом, в качестве числовой характеристики 
оптимальной оценки модели принимаем определитель матрицы С 


—1 
5ь =ае С =аеЦЕТЕ) (3.24) 
< 


Критерием выбора оптимального плана является минимизация 
чувствительности по всем возможным планам 


—1 
5 р =Шш5р = тп ае (ЕТ Е) (3.25) 
Зопт 5 
_ БТ 
т.к. определитель матрицы С = Е Ё всегда положителен. 


3.8 Синтез О-оптимальных тестирующих сигналов для 

идентификации динамических объектов 

Метод идентификации статистических объектов называется 
пассивная идентификация. 

Методы идентификации динамических объектов — активная 
идентификация — представляет собой сложную задачу, т.к. их 
применение связаны с необходимостью нарушения технологического 
режима идентифицирующего объекта и имеют значительно больше 
возможностей. 
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Активная идентификация предполагает подачу на вход объекта 
специального тестирующего сигнала При этом результаты 
идентификации, характеризуемые статистическими свойствами 
модели, полученной с помощью заданного метода оценивания 
(максимального правдоподобия, наименьших квадратов и т.п.), 
существенно зависят от структуры и параметров тестирующего 
сигнала. Выбор структуры параметров тестирующего сигнала 
представляют собой наиболее сложную часть общей задачи активной 
идентификации. 

Кроме структуры и параметров, включающих всю необходимую 
информацию о тестирующем сигнале, как функции времени, на 
качество идентификации модели влияет также часть общего времени 
идентификации, состоящая из временной предварительной пассивной 
идентификации возмущений и помех, действующих на объект, и 
временного активного воздействия на объект, определяемого длиной 
тестирующего сигнала. Статистические свойства модели зависят 
также от моментов регистрации выходной координаты объекта 
относительно начала активного эксперимента идентификации. 

Синтез (построение) р — оптимального тестирующего сигнала 
для линейных по параметрам динамических объектов сводится к 
следующему. 

Существует несколько видов наиболее часто используемых 
сигналов: ступенчатый, псевдослучайный (рис. 3.15 а,б), сигнал 
Дирека (рис. 3.16). 


Т=сопѕі 


а) б) 
Рисунок 3.15 — Тестирующие сигналы: а) ступенчатый тестирующий 
сигнал; б) псевдослучайный тестирующий сигнал 


Для идентификации динамических объектов применяют сигнал 
Дирека или 6(1) — дельта тестирующий сигнал (рис. 3.16). 
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Рисунок 3. 16 — Тестирующий сигнал Дирека ( 4(1)- дельта сигнал ) 


Подача 5() на входе динамической системы дает реакцию 
системы в виде так называемой весовой функции =(—т). Причем 
длительность переходного процесса 7 . 


Задача построения О — оптимального тестирующего сигнала 
для линейных по параметрам динамических объектов сводится к 
следующему: 

1) находят р — оптимальный план для соответствующей 
регрессионной модели; 

2) на основе полученного плана строят оптимальный 
тестирующий сигнал. 

Предположим, что объект является линейным или допускает 
линеаризацию в пределах малых отклонений от некоторого 
статистического состояния. Пусть выход У(г) измеряется дискретно с 


некоторым шагом ДЕ, а вектор шума е” = (е[0],е[А!|,е[2Аи...,е[пА!). 
Входной сигнал х0) является кусочно-постоянной на интервале Д{ 
функцией, тождественно равной нулю при отрицательном времени (до 
начала эксперимента). 

В качестве модели линейного динамического объекта может 
быть принята регрессионная модель, описывающая поведение объекта 
в отклонениях от некоторого статического состояния: 


1-1 
Ул! |= У [тл]. хп! — (т + ПАт|-+е[пАи|, (3.26) 
т=0 
где о[тлт| — дискретная весовая функция (набор ординат); 
хп — (т-+ 11| —щ  дискретное значение отклонения входного 


сигнала от их значений в исходном статическом состоянии; У [МАЕ — 
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дискретное значение отклонения выходного сигнала; лл] — 
случайная помеха; [Ат = Т, = х[иА!-(т+1)Ат| — априорная оценка 
времени переходного процесса объекта. 

Выражение (3.26) показывает, что единственным сигналом, 
которым можно варьировать в процессе эксперимента, являются 
дискретные значения входного сигнала. Эти сигналы и берутся в 
качестве факторов планирования эксперимента. 

Таким образом, при Ат = Лі факторами планирования будут 

х[(п –2)л:] Е=1.2,...1 (3.27) 

Тогда, при условии, что х(г)=0 при г<, пространство 

планирования Х, будет ограничено системой неравенств 
-1< (п -)л |1 (3.28) 

Далее, известно, что для моделей типа (3.26) ортогональные 
планы, построенные с помощью ПФЭ или ДФЭ, являются 
одновременно и Ш) — оптимальными. Поэтому для получения 
оптимального тестирующего сигнала проще использовать ПФЭ или 
ДФЭ. При этом следует иметь в виду, что в выражении (3.26) 
отсутствует свободный член. И эта особенность позволяет разбить 
матрицу ПФЭ на две, каждая из которых является р — оптимальным 
планом для данной модели. Разбиение матрицы на две делают с 
учетом того, что в каждой из полученных матриц один из факторов 
должен быть зафиксирован на одном из уровней +1 или -1. 

Учитывая, что всего факторов [, то существует / 
симметричных пар эквивалентных Д- оптимальных планов для 
модели (3.26), полученных на основе плана ПФЭ типа 2’. Есть 
возможность найти и другие р оптимальные планы. 

Покажем правило нахождения всех возможных р — 
оптимальных планов для модели (3.26), матрицы которых имеют 2"" 
строк. Берется матрица ПФЭ 2'. Строки этой матрицы формируются 
по две в каждой из №=2` групп. В каждую группу включаются 
только симметричные строки матрицы ПФЭ. Затем все возможные Р 
— оптимальные планы получаются путем включения одной из строк 
каждой группы в матрицу плана. Тогда общее количество возможных 
р — оптимальных планов для модели (3.26), полученных на основе 
ПФЭ типа 2' будет равно 2”. 

Пример. Пусть необходимо найти четыре ординаты весовой 


функции, т.е. [=4. Тогда ПФЭ будет иметь план 2'=2“ =16 (табл. 
3.13). 
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Таблица 3.13 


МПЭ ПФЭ 2* 





Перестраиваем матрицу 2*, так чтобы в каждую группу были 
включены симметричные строки (табл. 3.14). р — оптимальные 
планы 2' = 2° =8 получаются путем включения восьми строк, взятых 
по одной из каждой группы. Общее количество р — оптимальных 
планов будет равно 2” = 2° = 256, т.к. количество групп № =2*' =8. 


Таблица 3.14 
МПЭ ПФЭ 2* с симметричными строками 
Номер Номер 


опыта опыта 
табл. 3.14 табл.3.13 





ОБ обретет С ЗОВ Пе НЙ КО ИИ ИЕ СНЫ ЕЕ ТЕН 
РИ ее Е 

Теперь возникает задача построения самого р — оптимального 
тестирующего сигнала. Для этого существуют два пути. 

Первый, наиболее простой, состоит в стыковке строк матрицы 
ЮО — оптимального плана. Полученный таким путем тестирующий 
сигнал имеет максимальную длину, причем измерение выходной 
величины начинается через / — тактов, начиная с ( + 1). 

Второй путь состоит в ”сшивании" строк на у тактов. 
Очевидно, что минимальная длина тестирующего сигнала достигается 
в том случае, если все строки матрицы "сшиваются" на у = (1-1) 
тактов. В этом случае измерение выходной величины осуществляется 
на каждом такте, начиная с (/ +1). Такое "сшивание" можно считать 
корректным, если получающийся при этом план является Э — 
оптимальным. Для нашего случая полученный путем "сшивания" 
план должен быть представлен одним из 256 планов. 

Используем метод "сшивания" для получения ЮР — 
оптимального плана 2'' = 2*' на основе матрицы таблицы 3.14. 


На основе таблицы 3.15 получили минимальной длины ОЮ — 
оптимальный тестирующий сигнал +/+ /-/-/ -1-1+1-1 +1+1-1. 


Таблица 3.15 
— оптимальный план а 


Номер 
опыта 





В таблице 3.16 также представлен Рр-оптимальный план для 
нахождения четырех ординат весовой функции, однако длина 
тестирующего сигнала будет уже не оптимальной +/+1-1-1 +1-1-1-1 
-1+1-1+1. 
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Таблица 3.16 


0 — ТЕ план ДФЭ 


Номер Номер 
опыта опыта 





На рисунке 3.17 а, б) представлены сами тестирующие сигналы 
полученные соответственно по таблицам 3.15, 3.16. 


а) б) 
хит А и (1 



































| 
Ху № 


Рисунок 3.17 — Оптимальный тестирующий сигнал 


С учетом полученного оптимального тестирующего сигнала, 
модель линейного динамического объекта для каждого опыта примет 
следующий вид: 

У[АЕ |= охо + 1х1 +822 + 83Х 3; 
У! |= Вох + 810 + 82 + 83Х—2; 
УЗ! | = 0055 +810 + 8х0 + 3х1; 
Үүл] = 0003 +0 +290 +830; 

ү|5л|= 804 +8153 + 8х5 + 9321; 


Ү[6л]= 80%5 + 8х4 + 8223 + 8300; 
У А: |= вох + 81%5 + 8254 + 8323; 
УВ! |= 30%7 + 8156 + 8225 + 8354; 


Ул: = вохз + 8155 + 826 + 8325. 
Так как оптимальный тестирующий сигнал, построенный на 
основе таблицы 3.16, длиннее на один такт, то и реакций системы 
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будет получено на одну больше — не восемь, как для оптимального 
сигнала по таблице 3.15, а девять. 


Контрольные вопросы 

1. Каким образом при планировании эксперимента проводят 
выбор интервалов варьирования, выбор основного уровня плана 
эксперимента, выбор числа уровней? 

2. Для каких видов моделей проводят полный факторный 
эксперимент? В чем заключается принцип рандомизации плана 
эксперимента? 

3. Перечислите и сформулируйте свойства ПФЭ 2^. 

4. Для каких моделей возможно проведение дробного 
факторного эксперимента? 

5. Каким образом проводят планирование экспериментов при 
построении полной квадратичной модели? 

6. Дайте определение понятия «звездные точки». 

7. Какие критерии оптимальности планов Вы знаете? 

8. В каких случаях проводиться разбиение матрицы 2^ на блоки? 
Каким образом изменятся оценки параметров модели? 

9. Какие методы применяются для синтеза Ю — оптимальных 
тестирующих сигналов для идентификации динамических объектов? 
В чем они заключаются? 
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4. Экспериментирование 


4.1 Статистические методы исследования объектов 
4.1.1 Общие понятия и определения. Выборочная статистика 


Статистические методы используют идеи и методы 
многомерной математической статистики. Формализация процессов 
здесь осуществляется эмпирическими методами [11, 13]. 

Основной задачей математической статистики является 
разработка методов получения научно обоснованных выводов о 
массовых явлениях и процессах из данных наблюдений или 
экспериментов. 

В математической статистике вводится совсем новое и глубоко 
абстрактное понятие генеральной совокупности. Она представляет 
собой совокупность всех мыслимых наблюдений над случайной 
величиной при заданных условиях эксперимента. Если в результате 
эксперимента получены значения ху, х›...х,, то эти величины 
образуют выборку, а число п называется объёмом выборки. 

Термин выборочная статистика обозначает некоторое 
числовое значение, подсчитанное на основе обработки данных 
выборки. Следует помнить, что выборка охватывает не всю 
генеральную совокупность, поэтому выборочные статистики 
являются случайными величинами. Выборочными статистиками 
являются : среднее значение, выборочная дисперсия и выборочный 
коэффициент корреляции случайной переменной величины х. 


Р(х)= Е (х) — это плотность вероятности или плотность случайной 


х 
величины, где Ё (х) = | Та (х)ах — функция распределения. 
—00 
Выборочные статистики рассмотрим после функции 
распределения. 


4.1.2 Функция распределения выборок 


Функцией распределения Г (х) случайной величины & 
называют Е(х)= Р(& < х). 
Полученные при эксперименте числа ху, х,...х, считают 


совокупностью значений некоторой конечнозначной величины х, . 


При этом считают, что все эти числа являются различными 
элементами, независимо повторяются они или нет, т.е. каждый 
элемент выборки появляется лишь в результате одного измерения. 


Поэтому каждому элементу х, приписывают вероятность —. 
п 
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Полученное равномерное распределение величины х, называют 


эмпирическим или выборочным распределением. При увеличении 
объема выборки и распределение случайной величины х 


п 
оказывается близким к распределению генеральной совокупности х. 
Это утверждение доказывается основной теоремой математической 
статистики. 

Согласно этой теоремы: при и оо с вероятностью, равной 
единице, разность между функциями распределения случайных 


величин х, и х стремится к нулю 
т |Р, (0) Е(х) >0 

Практически при достаточно большом объеме выборки 
функцию распределения генеральной совокупности можно 
приближенно заменить выборочной функцией распределения. 

Вспомним, как вычисляется функция распределения. Для 
любого х, функция распределения равна сумме вероятностей 
значений величины х, не превосходящей х. 

Считают, что все элементы выборки имеют одинаковую 


1 
вероятность равную —. Поэтому функция распределения выборки в 
п 


каждой точке равна числу элементов выборки, меньших, чем х,, 


деленному на объем п. 
Пример. Пусть выборка состоит из элементов (-1.8; -1; 0.3; 2; 5). 
Объем выборки п=5. Тогда Ё (х) каждого элемента: 


ЕТ 
Е(х) = г = Е =0,2. Определим функцию при х=1 Е; (1) = =. при 


4 
х=4 В (4) = 5 = 0,8. Графическая иллюстрация примера приведена 


на рисунке 4.1. 


Е(х) 
1.0 
0.8} 


0.6 
0.4] 
А 


ЕЕ 10.3 _ = С 
-2 -1 0 1 2 3 4 5 х 








Рисунок 4.1 — Графическая иллюстрация примера 
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4.1.3 Числовые характеристики распределения 


Среднее выборки. Математическое ожидание 
Среднее выборки является наиболее эффективной выборочной 
статистикой. С помощью его оценивается математическое ожидание 
соответствующей случайной величины х, те х->мМ [х]. 
Математическое ожидание для бесконечного ряда совокупности 


00 
определяется следующим выражением: М [х|= [х./ (х)ах, где /(х)- 
—00 
плотность вероятности. Математическое ожидание для конечного 
п 
ряда совокупности – М [х|= Уху (х). 
к=] 


Случайная величина х, конечнозначна и имеет равномерное 


распределение. Поэтому среднее выборки есть среднее 
арифметическое элементов выборки 


х. 
х +Х,+...+Х, р> 





х = 
п п 

При достаточно большом п среднее может быть приближенно 
равным математическому ожиданию. При ограниченном числе 
опытов среднее выборки является случайной величиной, но все же 
оно дает более точное представление о математическом ожидании 
случайной величины х. 

При моделировании идеальное значение любой задаваемой 
величины при расчетах интерпретируется как математическое 
ожидание. 

Выборочная дисперсия 
Выборочная дисперсия случайной величины х, служит 


наилучшей оценкой дисперсии генеральной совокупности Х. 
Выборочная дисперсия также является случайной величиной и 
определяется выражением: 


п Е 5 
Б ба боа баа В 1-Х) 


п п 


Однако, учитывая, что на практике значение х не точно равно 
МЫ], те. х= М|х], дисперсия рү, ] дает, как правило, заниженную 





оценку рассеяния значений генеральной совокупности. Приведем пока 
без доказательства выражение: 
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га "1р. 


Из этого выражения видно, что истинный результат рассеяния 
значений генеральной совокупности [х] не совпадает с результатом 
рассеяния значений конечнозначной случайной величины Рх]. 


Значение р>, | получается несколько меньше значения Рх]. Поэтому 


п 


0х, | называют смещенной оценкой дисперсии Вх]. В качестве 


дисперсии выборки рассматривают величину 
п 


(а -х) 
В рУ х= х 
52 = рх или 52-1 НИИ (4.1) 
п 1 п 1 
Величина 5° является несмещенной оценкой генеральной дисперсии 
5? > рх] 
Докажем несмещенность выборочной дисперсии $7. 
Преобразуем выражение (2.1) 


(п—1)5? =У(х,-ху (4.2) 
Представим величину (х, – х) в следующем виде: 
(х; –х)= (х, -т,)-(х-т,) (4.3) 


здесь т, = МГ]. 
Подставим правую часть выражения (4.3) в (4.2): 


(и 05? = [м т) (ет) = 
= УС, -т, }* –2У`(х; -т, {х-т,)+ У(=-т,) = 
“Уст, зяби, Риби, = 


= (к-т, -п(х-т,)°. 


Найдем математическое ожидание обеих частей выражения 


(4.4): 


(4.4) 





п 


мст д = или 


= пр[х]- „Рі = Вх п –1) 


(„-1)м [52 | („ рр]. 
Отсюда следует, что математическое ожидание выборочной 
дисперсии 5” равно генеральной дисперсии др]. Значит 8° является 
несмещенной оценкой генеральной дисперсии БЫ]. 
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При моделировании идеальное значение любой задаваемой 
величины середнеквадратического отклонения при расчетах 
интерпретируется как дисперсия. 

Выборочный коэффициент корреляции. 
Понятие о корреляции 

В математическом анализе зависимость между двумя 
величинами выражается функцией. Эта зависимость однозначно 
связывает независимую переменную с другой зависимой переменной. 
Эту зависимость называют функциональной. Ее устанавливают с 
помощью строгих логических доказательств на основе определенных 
физических законов. 

Намного сложнее дело обстоит с понятием зависимости между 
случайными величинами. Например, если при изменении случайной 
величины х изменится и другая случайная величина у, то нельзя 


сразу ответить, является ли это изменение результатом зависимости 
от х, или оно обязано лишь влиянию случайных факторов. 

В природе существует причинно-следственная связь между 
случайными величинами. Эта связь особого рода, при которой с 
изменением одной величины меняется распределение другой. Такую 
связь называют стохастической связью. Оценка силы этой связи 
представляет собой сложную задачу математической статистики. В 
общем виде эта задача до сих пор не решена. Однако существуют 
показатели, оценивающие те или иные стороны стохастической связи. 
При этом основным показателем является коэффициент корреляции. 

Необходимым и достаточным условием корреляции между 
случайными величинами служит неравенство 


Му -т, х= т, )1=0, 


где тт —— математическое ожидание величины х и у 


у 
соответственно. 


2 
Величина М [С —т, Хх И, у = бух Носит название 


корреляционного момента или чаще называют ковариацией. 
На практике удобнее пользоваться безразмерным 
коэффициентом 
2 
З м у Хх т.) = бух 


О Ы 
| Рур; 68; 
Этот коэффициент называется коэффициентом корреляции. 
Для независимых величин коэффициент корреляции равен нулю. 
Коэффициент корреляции характеризует не только наличие, но 
и силу стохастической связи между уи х. Максимальная корреляция 





равна р= +1. 


ЭЗ 


При р>0 случайные величины у и х одновременно возрастают 
или убывают, а при р<0 с возрастанием одной величины другая 


убывает или наоборот. 

Коэффициент корреляции, как показатель зависимости 
случайных величин, имеет и слабую сторону: из равенства о = 0 еще 
не следует независимость величин у и х. Величины могут быть даже 
связаны функционально, а коэффициент корреляции может быть 
равен нулю. Это возможно при нелинейной зависимости. Однако, в 
тех случаях, когда заведомо известно, что между величинами 
существует линейная связь коэффициент корреляции является 
достаточно полным показателем. 

Выборочный коэффициент корреляции г является лучшей 
оценкой коэффициента корреляции р (генерального коэффициента 
корреляции). 

Пусть проведено т испытаний и при каждом из них 
фиксировались значения двух случайных величин. Тогда получается 
т- пар выборочных значений (х, у, ); (х, уз)...(х,„,у„): Эти точки 


можно рассматривать на плоскости и их расположение дает нам 
представление о силе корреляции. 





а) 0) 8) 


Рисунок 4.2 — Расположение экспериментальных точек на плоскости 


На рис. 4.2 приведены примеры, соответствующие сильной 
корреляции, случай а); слабой 6) и полному ее отсутствию, случай в). 

Выборочный коэффициент корреляции находится аналогично, 
как и генеральный коэффициент корреляции. 

Выборочный корреляционный момент равен 

2 | Е Е 
Ра хее) 
т-1:= 
а выборочный коэффициент корреляции: 
Рух 


=, 
8285 
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й 
где 5 - 9. в выборочные дисперсии 


1 2 2 Б е 
о 52 = —— (>, =). 


т=1 т- 1 т=1 
Выборочный коэффициент корреляции с увеличением объема 
выборки 7 приближается к коэффициенту корреляции р. 


Нужно помнить, что выборочный коэффициент корреляции дает 
оценку для коэффициента корреляции р только в тех случаях, когда 
совместное распределение двух случайных величин уи х является 
нормальным. При этом зависимость между у и х будет линейной. 
Однако, если только каждая случайная величина имеет нормальный 
закон распределения, то эти величины могут быть 
некоррелированными, но зависимы друг от друга. Только совместное 
нормальное распределение случайных величин дает возможность 
использовать выборочный коэффициент корреляции г для анализа 
стохастической связи между величинами. 


4.1.3 Квантили распределения 


Рассмотренные выборочные параметры: среднее, среднее 
квадратическое отклонение и коэффициент корреляции выборки 
являются приближенными оценками соответствующих генеральных 
параметров. Погрешность этих оценок будет тем меньше, чем больше 
объем выборки. Есть способы, с помощью которых можно оценить 
саму погрешность. Для этого переходят от точечных оценок 
параметров к оцениванию доверительных интервалов параметров. 
При получении интервальных оценок часто используют так 
называемые квантили. 

Квантилем, отвечающий заданному уровню вероятности Р, 


называют такое значение х=х, , при котором функция 


распределения Ё (х) принимает значение, равное Р, т.е. Е [х,]= Р, 
где Р — заданный уровень вероятности. 


Другими словами квантиль есть такое значение случайной 
величины х, при котором Е[х<х,]= Р. 


Вероятность Р, задаваемая в процентах, дает название 
соответствующему квантилю, например ху, называется 40%-ым 


квантилем. 
Квантили стандартного нормального распределения 
(распределение с параметрами т, = 0, с =1) обозначаются буквой и,. 


1 
Они сведены в соответствующих таблицах. Если р < Э? то подбирая 


ЭЭ 


1 1 
5 х)= = - и =—Х. = 
такое х, для которого Ф(х) 5 р и находим и, Если р> - 


1 
то подбираем такое х, для которого Ф(х)= РЭ и тогда и, = х. 


Например, 40 % квантиль будет равен из. = -0.25, 85 %, квантиль 
Квантиль общего нормального распределения о, с параметрами 


т, и с выражается через квантиль и,: 
ор = т. +Си. (4.5) 


Если известны два квантиля случайной величины х, ИХ, ТО 
Ех =х<х,]=9- Рр, 
Понятие квантиля используется не только для нормального, но и 


для большинства встречающихся распределений. 


Квантиль х у называется медианой распределения. Если 
2 


распределение случайной величины симметрично, то х и = М (х). 


Например, распределение случайных ошибок симметрично. 
Поэтому для этого распределения можно использовать как 
математическое ожидание, так и медиану. 


4.1.4 Получение интервальных оценок 


Все выборочные параметры являются случайными величинами, 
а, следовательно, и их отклонения от генеральных параметров, т.е. 
погрешности, также будут случайными. Поэтому оценка этих 
погрешностей носит вероятностный характер. Можно лишь указать 
вероятность той или иной погрешности. Чтобы решить подобную 
задачу нужно найти вероятность того, что отклонения выборочного 
параметра у от исследуемого генерального Ө не превосходит по 
абсолютной величине некоторого заданного числа є, т.е. находятся в 
пределах от —= до +=. Обозначим отклонение через Ди. 
Поставленную задачу можно легко решить, если известна функция 
распределения Е(х) или плотность распределения у(х) величины Лу 


РЇАи< 2] Е(6) (е) | (а (4.6) 


Распределение Лу иногда удается точно определить по 
элементам выборки. В других случаях это распределение зависит 
только от объема выборки п и его можно вывести теоретически. Если 
бы при этом было известно математическое ожидание выборочного 
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параметра у, то разность 
М["]- МГАу|= М[У]- Му 60]= Му] М[У]+0=0 
дает точное значение генерального параметра. К сожалению М [у], 
как правило, не известно. 
Определим генеральный параметр. Находят по выборке одно 
значение у, выборочного параметра у и принимают его за 


приближенное значение генерального параметра 0. Затем, используя 
выражение (4.6) оценивают это приближение. Действительно, 
задаваясь некоторым положительным числом =, находим вероятность 


того, что|Ли| =|-6 <=. Но так как Уд есть одно из допустимых 


значений выборочного параметра у, то вероятность неравенства 
|ио — 6] < = также равна этой вероятности. Отсюда получаем формулу 


Ри, -9]- в(е)-Е(е)= | бах 


которая позволяет сравнивать найденное значение выборочного 
параметра с неизвестным генеральным параметром. 

Неизвестный генеральный параметр можно представить в виде 
ио -= << +. 

Это неравенство отличается тем, что неизвестная величина Ө, 
которая является неслучайной величиной, оценивается случайными 
границами, т.к. выборочный параметр является случайным. 

Таким образом, любая статистическая оценка есть оценка вида 


у <б<и, 
где и ‚и некоторые случайные величины. 


Придавая у ‚у конкретные значения, вычисляем вероятность 
соответствующей оценки. 

В качестве границ у,’ наиболее удобно брать квантили 
случайной величины У. При обработке наблюдений для оценок 
генерального параметра берут симметричные квантили. В этом случае 
вероятности Р соответствует оценка уу, <0<у1,, 

5 2 

Рассмотрим, как пользоваться оценками на практике. Для этого 
дадим следующие определения. 

Событие называется абсолютно достоверным, если оно 
появляется при любом осуществлении комплекса основных факторов. 
Абсолютную достоверность нельзя установить никакой самой 
длительной проверкой. Ее можно вывести лишь теоретически, путем 
логических умозаключений. Сюда относят обычно математические 
истины. 

Большинство же привычных достоверных событий при 


Эл 


рассмотрении не является абсолютно достоверным. Например, нельзя 
считать абсолютно достоверным тот факт, что подброшенная монета 
упадет гербом или числом, т.к. у монеты есть и другие состояния 
равновесия - это ребро. Однако, в данном случае можем утверждать, 
что монета упадет либо гербом, либо числом. Такая достоверность 
называется практической достоверностью. 

Использования принципа практической достоверности 
позволяет не доводить вероятность оценки до единицы. Принимаемый 
при этом уровень вероятности называется доверительной 
вероятностью. В зависимости от конкретных обстоятельств в 
качестве доверительной вероятности берут обычно значения: 0.95 = 
0.99; реже 0.90;0.999. 

Соответствующие доверительной вероятности квантильные 
границы называются доверительными границами, а образуемый ими 
интервал — доверительным интервалом или еще называют 
доверительной оценкой или интервальной оценкой. 

Величина @=1—р называется уровнем значимости. Уровень 


значимости соответствует практически невозможному событию. 
Уровень значимости и уровень достоверности в сумме дают единицу. 
Обычно значения уровней значимости берутся в пределах: 0.05 +0.01 
и реже 0.1; 0.001. 


4.2 Генераторы случайных чисел. Псевдослучайные числа 

При построении стохастических имитационных моделей 
необходимо обеспечить возможность учета стохастических 
воздействий. Результаты статистического моделирования 
существенно зависят от качества выбранных исходных 
последовательностей случайных чисел. На практике используется три 
основных способа генерации случайных чисел: аппаратный 


(физический), табличный (файловый) и алгоритмический 
(программный) [6, 13]. 
Определим основные требования, которым должен 


удовлетворять идеальный генератор. То есть псевдослучайные 
последовательности чисел должны: 


1) состоять из квазиравномерно распределенных чисел; 
2) содержать статистически независимые числа; 
3) быть воспроизводимыми; 


4) программа генерации должна иметь наименьшие затраты 
машинного времени и занимать минимальный объем памяти. 
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4.2.1 Методы получения псевдослучайных чисел 


Метод средних квадратов 

Одной из первых алгоритмических процедур получения 
псевдослучайных чисел был метод средних квадратов, который 
заключается в следующем. Например, для восьмиразрядного числа 
(половина равна 4 разрядам) берется начальное значение: 
хо = 0.2152, возводится в квадрат: 

(х)= 0.04631104, берется “среднее” число: 
хт = 0.6311 и снова возводится в квадрат: 

(х1 )? = 0.39828721, 

х = 0.8287 ит. д. 

Однако этот метод работает неудовлетворительно: существует 
наличие корреляции между числами, а иногда и отсутствие характера 
случайности. 

Конгруэнтные методы 

Самое широкое применение при моделировании на ЭВМ 
получили конгруэнтные методы генерации псевдослучайных 
последовательностей, в основе которых лежит фундаментальное 
понятие конгруэнтности. 

Два целых числа @ и В конгруэнтны (сравнимы) по модулю т, 
где т — целое число, тогда и только тогда, когда существует такое 
целое число К, что 

а- В= Кт, 
то есть разность а - В делится на т, и если числа а и В дают 
одинаковые остатки деления на абсолютную величину числа т. 
Например: 
1984 = 4 ( п10410),5008 = 8 ( той 10°) ит.д. 

Величина т берется равной длине машинного слова т = О. гдер 
— число бит в машинном слове. 

Конгруэнтный метод представляет собой арифметическую 
процедуру для генерирования конечной последовательности 
равномерно распределенных чисел. С использованием нескольких 
рекуррентных формул было построено множество конгруэнтных 
алгоритмов. Среди них наиболее известными являются 
мультипликативный, смешанный и аддитивный методы. 

Мультипликативный конгруэнтный метод 

Основная формула этого метода имеет вид: 

хн = ах, (той т), 

где а и т — неотрицательные целые числа. 

Согласно этому выражению, должны взять последнее случайное 
число х, умножить его на постоянный коэффициент а и взять модуль 
полученного числа по т, то есть разделить на ах, и остаток считать 
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как х„1. Поэтому для генерирования последовательности чисел х; 
необходимо начальное значение хо, множитель а и модуль т. 

Любой генератор псевдослучайных чисел может дать лишь конечное 
множество чисел, после которого последовательность будет 
повторяться. Период, или длина последовательности Ё, зависит от 
ЭВМ и от выбранного модуля, а статистические свойства зависят от 
выбранного начального значения и множителя. Таким образом, 
выбирать а, хо и т следует так, чтобы обеспечить максимальный 
период и минимальную корреляцию между генерируемыми числами. 
Правильный выбор модуля не зависит от системы счисления. Обычно 
берут т равным длине машинного слова: 

— для двоичной ЭВМ т=2? ‚ где р - число двоичных цифр (бит) 
в машинном слове; 

— для десятичной ЭВМ т=/1 0“ ‚ где 4 - число десятичных цифр в 
машинном слове. 

Тогда, если правильно выбрать а и хо, то максимальный период 
будет равен 

[=2”?=пиИ4 когда 6>2 (двоичная система счисления); 
[=5 .10“?=т/20 когда 4>2 (десятичная система счисления). 

Для двоичной системы а = 8Т + 3, где Т может быть любым 
целым положительным числом, хо — любое положительное, но 
нечетное число. 

Для десятичной системы а = 200Т + О, где Т любое целое 
положительное число, О может принимать одно из следующих 
значений + (3, 11, 13, 19, 21, 27, 29, 37, 53, 59, 61); хо — любое 
положительное число, которое не делится на 2 и на 5. 

Процедура вычисления случайных чисел между 0 и 1: 

1) выбрать в качестве хо произвольное нечетное число; 

2) вычислить коэффициент а; 

3) найти произведение ахо, содержащее не более 2р значащих 
разрядов; 

4) взять Б(а) младших разрядов в качестве первого члена 
последовательности х1, а остальные отбросить; 

5) определить дробь ху = х; /2” (х =х /10“ ) из интервала (0, 1); 

6) приравнять хо = х1; 

7) вернуться к пункту 3. 

Пример. Надо получить последовательность чисел для случая 
Ь=4. 

1. Выбираем хо=0111(7). Берем Т=/. 

2. а=8Т-3=8-3=5, а=0101. 

3. ахо=(0101)(0111)=00100011. 

4. х 1 =0011. 
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5. Проверяем дробь интервала (0,1) х;=х/2*=3/16=0,1875. В 
результате получаем: 

1) ахо=(0101)(0111)=00100011, х,=0011, х=3/16=0,1875; 

2) ахт=(0101)(0011)=00001111, х›=1111, х=15/16=0,9375; 

3) ахг=(0101)(1111)=01001011, хз=1011, хз=11/16=0,6875; 

4) ахз=(0101)(1011)=00110111, х;=0111, х=7/16=0,4375. 

Таким образом, генерируемая последовательность чисел х, і = 0, 
1, ..., №, представляет собой равномерно распределенные на интервале 
(0, 1) случайные числа. 

Для дискретной последовательности используют 2" случайных 
чисел на интервале (0,1). Поэтому закон распределения такой 
дискретной последовательности называют квазиравномерным 
распределением. При этом случайная величина принимает значения: 


і оп 





1 1 . 
с вероятностью р, а. гдеі = 0, 1,..., 2" -1. 
1 


Тогда математическое ожидание и дисперсия дискретной 
последовательности соответственно равны: 


2" ; г 1 ры ЕО" 1 
> ЕТ 1) = "1127.2 в. 

















25 ЕО ИК 
0 Б) с а 12: атъ 
Для равномерной случайной последовательности этого же 
интервала (0,1) имеем: 


| 
м(&)= | х7 ()ах => ‚где /(х)- функция плотности, 
0 
== 
Јо) = | 


- мерла =>. 


При достаточно больших и свойства законов распределения 


квазиравномерной последовательности и равномерной 
последовательности совпадают. Кроме того, для получения значений х 
используется алгоритмический способ. Поэтому такие 


последовательности по своей сути являются детерминированными и 
называются исевдослучайными. 
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Смешанный метод 

Работа этих генераторов основана на использовании формулы: 

хн = (ах, + С) (той т), С>0. 

С вычислительной ТОЧКИ зрения смешанный метод генерации 
последовательности неотрицательных целых чисел сложнее 
мультипликативного на одну операцию сложения, но при этом 
возможность выбора дополнительного параметра позволяет 
уменьшить возможную корреляцию получаемых чисел. Однако 
экспериментальная проверка качества генерируемой 
последовательности чисел на основе этой формулы является сложнее. 


4.2.2 Проверка качества квазиравномерной последовательности 
псевдослучайных чисел 

Эффективность моделирования на ЭВМ существенно зависит от 
качества исходных последовательностей псевдослучайных чисел. 
Поэтому, прежде чем приступить к реализации моделирующих 
алгоритмов на ЭВМ, необходимо выполнить тестирование случайных 
чисел, которое включает в себя проверку равномерности, проверку 
стохастичности и проверку независимости. 

Проверка равномерности выполняется по гистограмме. Для 
проверки выдвигается гипотеза о равномерности распределения чисел 
в интервале (0, 1). Затем этот интервал разбивается на т равных 
частей. Это означает, что при генерации последовательности { х; } 
каждое из чисел х; с вероятностью 

Р,==.ј=1л 
попадет в один из частей интервала. Всего в каждый /-ый подинтервал 
попадет №; чисел, причем общее количество чисел: 
к 
мМ= УМ. 
Е 

Относительная частота попадания случайных чисел в каждый из 
подинтервалов равна №/ №. На рисунке 4.3 показан пример 
гистограммы распределения чисел. 

Пунктирная линия соответствует теоретическому значению Р}, а 
сплошная — экспериментальному №;/ №. Очевидно, что при большом 
М№ экспериментальная гистограмма приблизится к теоретической. 
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Рисунок 4.3 — Гистограмма распределения чисел 


Оценка степени приближения может быть проведена с 
использованием критериев согласия (Колмогорова, Пирсона, 
Стьюдента, Фишера и т.п.). На практике обычно принимается К = 
20 = 50, №= (10° = 10°) К. 

Проверка стохастичности 

Среди наиболее часто используемых тестов являются 
интервальный тест и автокорреляционные тесты. 

Интервальный тест проводится методами комбинаций и серий. 
Сущность метода комбинаций сводится к определению закона 
распределения длин участков между единицами (нулями) или закона 
распределения появления числа единиц (нулей) в п-разрядном 
двоичном числе х; На практике длину последовательности № берут 
достаточно большой и проверяют все и разрядов или только / старших 
разрядов числа х;. 

Теоретически закон появления ] единиц в / разрядах двоичного 
числа х; описывается биноминальным законом распределения: 

Р(}1) = СІР 1), 
где Р (7,1 ) — вероятность появления ] единиц в / разрядах числа х;, 
Р(1) = =Р (0) = 0.5 — вероятность появления единицы (нуля) в любом 
разряде числа х; сі — количество комбинаций появления единиц 
(нулей) в / разрядах: 
п 
(1-1! 

Тогда, если выбрать длину выборки № фиксированной, то 
теоретически ожидаемое число появления случайных чисел х; с ] 
единицами в проверяемых / разрядах будет равна 

П=МСИР' О. 

После нахождения теоретических и экспериментальных 

вероятностей Р(/,/) или чисел 7; при различных значениях / <и 


С 
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гипотеза о стохастичности проверяется с использованием критерия 
согласия. 
Проверка независимости 

Проверка независимости элементов последовательности {х,} 
проводится на основе вычисления корреляционного момента. 
Случайные числа & и у называются независимыми, если закон 
распределения каждого из них не зависит от того, какое значение 
приняло другое. Для независимых величин корреляционный момент 


К „= хх, -м(2)у, -м(@)р,. 
г) 
и коэффициент корреляции равны нулю 


= К 
Ра | 


где Р; — вероятность того, что (&,77) примет значение (х;, у;); 
С х,су — среднеквадратические отклонения величин ёи 7.. 


Таким образом, независимость элементов последовательности 
{х;} может быть проведена путем введения в рассмотрение 


последовательности {У;}, для которых у; =х; ‚ Где т — величина 


і+т 
сдвига последовательностей. 

При проведении оценок корреляции на ЭВМ удобно для 
вычисления использовать следующие выражения: 


1 > 2155 о 
х.х. — х 
РЕН №М-т 1=1 а (М-— = 1=1 1= 


Рас (р(х )р(х;,„) 


1 М№-т 2 1 №М-т 
Вх = 2х; { р 2 х; 
1=] Мт 1=1 


Б) У 
Вх; ТАРР > НЕ" Е са 














2 


1= 
При любом т * 0 для достаточно больших № с доверительной 
вероятностью р справедливо соотношение: 


Ра) 58 


где обычно принимает значение в интервале д=0,9-0, 99. 





Если найденное эмпирическое значение р 8, (г) находится в 
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указанных пределах, то с вероятностью Д можно утверждать, что 
полученная последовательность чисел х; удовлетворяет гипотезе 
корреляционной независимости. 


4.3 Моделирование случайных сигналов и процессов 
4.3.1 Формирование возможных значений случайных величин с 
заданным законом распределения 
Исходным материалом служат случайные величины с 
равномерным распределением в интервале (0, 1): 


1 0<х<1 
ал 6 ‚где /(х) — функция распределения. 


Есть простой метод преобразования равномерного 
распределения в требуемое. Его называют прямым методом. Он 
заключается в преобразовании случайных величин. Пусть переменная 
у является функцией от х: 

у= ф(х). 

Тогда для каждого значения х однозначно определяется у. 

Пусть х — случайная величина, тогда у будет так же случайной 
величиной. И если х имеет функцию распределения } (х), то у будет 
иметь так же функцию распределения } (у), но зависящую уже от 
(х) и ф(х). 

Известно, что для взаимно однозначного соответствия между х и 
у имеет место следующее равенство: 

Л б) ах=ђ 0) ду, (3.1) (4.7) 

то есть вероятность того, что случайная величина х находится в 
интервале между х; и х+ах, равна вероятности того, что 
преобразованная случайная величина у будет находиться в интервале 
у; и у+ау. 

Последнее выражение можно изобразить графически (рис. 4.4). 
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Рисунок 4.4 — Прямой метод преобразования 


Согласно выражению (4.7) находится требуемое распределение 
случайной величины х: 


0) ОБ: 


ф'(х 
ах 


Требуемое распределение № (у) является функцией от 
распределения случайной величины х и производной от функции ф 
(х). Абсолютное значение производной принято потому, что функция 
распределения не может принимать отрицательных значений. 
Функцию у = ф (х) можно найти из (3.1): 

П а Л ‚(х) .. Л „(х ) 
Ф\х)= или ф(х)= | 
Ў.) 200) 

Пусть р (у) = сопѕі. Тогда имеем равномерное распределение, и 

выражение (4.8) принимает вид: 


о 1 У (в или 96) = Е (х). 


Функция ф (х) для данного случая совпадает с интегральной 
функцией распределения Е; (х). Это выражение можно использовать 
для преобразования случайных чисел, то есть для реализации 








Ах (4.8) 
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некоторой операции над числом у, имеющим равномерный закон 
распределения, и преобразование его в число х, имеющее заданный 
закон распределения: 


0с зБ уб). 


Эта формула является рабочей. Она может быть использована 
для практических целей. 

Пример. Необходимо получить случайные числа х; с 
показательным законом распределения: 


Р(х)= 4.7 (х> 0). (4.9) 


Согласно этому выражению имеем: 


х; 
их 
у;=4 еза; 
—00 
Здесь у; — случайная величина с равномерным законом 


распределения в интервале (0, 1). 
После вычисления интеграла получим: 


4х; или х; = (У | 


о і 
И далее, учитывая, что случайная величина у; = 1-у имеет так же 
случайное распределение в интервале (0, 1), запишем 
1 
ОЕ (4.10) 


1 


Однако для многих практических случаев, например для 
нормального распределения, уравнение (4.9) не имеет точного 
решения. Кроме того, даже простое на первый взгляд выражение 
оказывается неудобным для решения его на ЭВМ. Поэтому на 
практике обычно пользуются приближенными приемами 
преобразования случайных чисел, основанными на кусочной 
аппроксимации функции плотности (рис.4.5). 


4.3.2 Алгоритм формирование случайных величин 

Пусть требуется получить последовательность случайных чисел 
{ у; } с функцией плотности ў, (у), возможные значения которой лежат 
в интервале (а,в). Представим }, (у) в виде кусочно-постоянной 
функции. Для этого разобьем интервал (а, в) на т интервалов и будем 
считать, что ў, (у) на каждом интервале постоянна. Чтобы 
аппроксимировать /, (у) наиболее удобным способом для 
практических целей, целесообразно разбить (а, в) на интервалы так, 
чтобы вероятность попадания случайной величины 7 в любой 
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интервал (ак, аці) была постоянной, то есть не зависела от номера 
интервала К. Тогда для вычисления а; можно использовать выражение 


ак] 1 
е (4.11) 
Ч 


Случайную величину 7 можно представить в виде 7 = ах + Пь 
где ак — абсцисса левой границы К-го интервала; 7. — случайная 
величина, возможные значения которой располагаются равномерно 
внутри К-го интервала. 


НЫ 
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р” 






































а а 


К К+1 


Рисунок 4.5 — Кусочная аппроксимация функции плотности 


Алгоритм этого способа получения случайных чисел сводится к 
следующим действиям: 

1) генерируется случайное равномерно распределенное число х; 
из интервала (0, 1); 

2) с помощью этого числа случайным образом выбирается 
интервал (ах, ак); 

3) генерируется число хь: и масштабируется с целью 
приведения его к интервалу (а, акі), то есть помножается на 
коэффициент (аы - ак) хыт; 

4) вычисляется случайное число у; = ак + (ан - а) ха с 
требуемым законом распределения. 

Для выбора интервала (ак, ар) с помощью случайного числа х; 
строят таблицу (формируют массив), в которую предварительно 
помещают номера интервалов К и значения коэффициента 
масштабирования, определенные из соотношения (4.11) для 
приведения числа к интервалу (а, в). Получив случайное число х;, с 
помощью этой таблицы определяют абсциссу левой границы а; и 
коэффициент масштабирования (а, - а. 

Достоинство этого приближенного способа преобразования 
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случайных чисел состоит в том, что при реализации на ЭВМ 
требуется небольшое количество операций для получения каждого 
случайного числа, так как операция масштабирования выполняется 
только один раз перед моделированием, и количество операций не 
зависит от точности аппроксимации, то есть количества интервалов т. 

Можно привести еще способы преобразования равномерно 
распределенных случайных чисел { х; } в последовательность с 
заданным законом распределения { у; } на основе предельной теоремы 
теории вероятностей. 


4.4 Центральная предельная теорема 

Если &, &, .., бү — независимые одинаково распределенные 
случайные величины, имеющие математическое ожидание а и 
дисперсию о”, то при № -> оо закон распределения суммы случайных 

М 
чисел У; х; неограниченно приближается к нормальному закону. 

1=1 

Рассмотрим применение предельной теоремы на примере. 

Пусть требуется получить последовательность случайных чисел 
{ у }, имеющих нормальное распределение с математическим 
ожиданием а и средним квадратичным отклонением о. 


(у-а} 


е 20? 
2ло 


Случайные числа у; формируются в виде сумм 
последовательностей случайных чисел { х; }, имеющих равномерное 
распределение на интервале (0, 1). Тогда можно использовать 
центральную теорему. 

Если независимые одинаково распределенные случайные 
величины х], х2, ..., ху имеют каждая математическое ожидание а! и 

М 
среднее квадратичное отклонение с, то сумма Ух; асимптотически 
1=1 
нормальна с математическим ожиданием а = № и с=с Ум. 
М 
Расчеты показывают, что сумма У, х; имеет распределение, близкое к 
1=1 
нормальному, уже при № = 8 - 12. 
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4.5 Правило автоматической остановки имитационного 

эксперимента 

Первый способ. Простейший способ состоит в задании 
требуемого количества реализаций №. Однако такой подход является 
грубым из-за того, что на этапе тактического планирования 
неизвестны распределения выходных переменных. 

Второй способ состоит в задании доверительных интервалов 
для выходных переменных. Остановка прогона машинной модели 
происходит при достижении заданного доверительного интервала. 

Третий способ выполняется путем двухэтапного проведения 
прогона, когда сначала делается пробный прогон из №; реализаций, 
позволяющий оценить необходимое количество реализаций №. Если 
№> М , то прогон можно закончить, в противном случае необходимо 
набрать еще (/- №) реализаций. 

Четвертый способ выполняется путем последовательного 
анализа для определения минимально необходимого количества 
реализаций №, которое рассматривается при этом как случайная 
величина, зависящая от №-/ предыдущих реализаций при достижении 
нормального закона распределения. Такой подход позволяет 
проводить анализ результатов эксперимента на базе характеристик 
выборок. 

Рассмотрим четвертый способ на примере двухфакторной 
линейной (регрессионной) модели (теоретически данный вопрос 
рассматривался в п.п. 5.2.2). 


зА 











х 


Рисунок 4.6 — График эксперимента 


На рисунке 4.6 показаны точки (х,у;), 1=1, №, полученные в 
машинном эксперименте с моделью системы. Делаем предположение, 
что модель результатов машинного эксперимента графически может 


быть представлена в виде прямой линии У=БЬ +Бх, где У - 
величина, предсказываемая регрессионной моделью. Требуется 
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получить такое значение коэффициентов В, и В, при которых сумма 


квадратов ошибок является минимальной. 
На рисунке 4.6 ошибка е, для каждой экспериментальной точки 


определяется как расстояние по вертикали от этой точки ДО ЛИНИИ 


регрессии У. 


^ 


е; =У, —У, =Б+Ьх, -, 


І 


1 1 


№ 
Ф= УВ, + рух; — У, ) — сумма квадратов ошибок, 


ОФ М М 
дро 1=1 1=1 
ОФ 


м и м 
Ору і=1 і=1 і=1 
Решая систему нормальных уравнений получим 


Ь УУ -Ужу хи, 
0 ТЕ э 2 , 
мух -| Ух 
і=1 і=1 
= М, х1, БОО 
у м м \2’ 
му |524 
= = 


№ — число реализаций при моделировании системы. 














Мерой ошибки регрессионной модели служит 
среднеквадратическое отклонение 
М Я М А 
26. У(6,+6,Х,-Ү,) 
с, = і=1 == і=1 Я 
№М-2 М-2 


Для нормального распределенных процессов приблизительно 
67% точек находится в пределах одного отклонения о, от линии 
е 


регрессии — труба А, и 95% точек — в пределах 20`,- труба В (рис. 
4.7). 
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Х 


Рисунок 4.7 — Пределы отклонения экспериментальных точек 


Для проверки точности оценок Р.,Ь и адекватности модели 
используются критерии Стьюдента и Фишера. 


Контрольные вопросы 

1. Перечислите статистики генеральной совокупности и 
характеристики выборки. 

2. Дайте определение функции распределения и постройте 
данную функцию для выборки. 

3. Каким образом при экспериментировании соотносятся 
среднее выборки и математическое ожидание? 

4. Каким образом при экспериментировании соотносятся 
среднее квадратическое отклонение и дисперсия? 

5. Дайте определение корреляции, корреляционного момента и 
коэффициента корреляции. 

6. Объясните понятие квантили распределения. 

7. Объясните принцип получение интервальных оценок. 

8. Каким образом получают оценку математического 
ожидания? 

9. Каким образом получают оценку генеральной дисперсии? 

10. Перечислите известные Вам методы получения 
псевдослучайных чисел. Подробно раскройте метод средних 
квадратов. 

11. Подробно расскажите о получении последовательности 
случайных чисел мультипликативным конгруэнтным методом и 
смешанным методом. 

12. В чем заключается проверка равномерности 
последовательности случайных чисел? 
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13. В чем заключается проверка независимости 
последовательности случайных чисел? 

14. В чем заключается проверка стохастичности 
последовательности случайных чисел? 

15. Каким образом проводят формирование возможных 
значений случайных величин с заданным законом распределения? 

16. В чем заключается алгоритм получения случайных величин? 

17. Какие правила автоматической остановки имитационного 
эксперимента Вы знаете? 
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5. Анализ результатов 


5.1 Оценка результатов наблюдений 
Целью статистической оценки является отыскивание оценки 


генерального параметра — идеального значения на основе 
выборочного параметра — экспериментального значения. Будем 
предлагать, что наблюдаемая случайная величина имеет нормальное 
распределение. 


Основным оцениваемым параметром является генеральное 
среднее, т.е. математическое ожидание. 
5.1.1 Оценка математического ожидания 
Математическое ожидание легко определяется, если известна 


генеральная дисперсия С°. Зная с’, можно дать оценку для 
математического ожидания даже по одному наблюдению. 

Пример 1. Определить доверительный интервал 
математического ожидания "а" по одному наблюдению х=3, если 
известно, что генеральная совокупность имеет нормальное 
распределение с значением дисперсии с=0.9. В качестве 
доверительной вероятности возьмем Р=0.96. 

Берем в качестве доверительных границ симметричные 
квантили 
<а<и, или иу 0,96 =9 0110096, 

2 Е) Е абы. 
Тогда имеем: иу < а < идоз; по таблице – 2.06 < а < 2.06. 
Для общего нормального распределения из (4.5) 
З+ о -ио02 <а<З+ о: ио оз 


3+0.9(-2.06)<а<3+0.9.2.06 


1.146 <а<4.854 
Если над случайной величиной хпроведено несколько 
наблюдений, то для оценки математического ожидания можно 


использовать выборочное среднее х и тогда 
Х+05 ‘Ир <а<х+о; "Шр 
2 2 


о 
Если генеральная дисперсия известна, то со; Е. и 
п 


22 ЕЕ => 
окончательно получаем: Х + =: <а<х+-=‘Инр, 
Ча “ЕР д 
2 2 
Отсюда видно, что уменьшение интервальной оценки обратно 
пропорционально корню квадратному из числа наблюдений. 
т! 


Пример 2. Оценить математическое ожидание "а" по 
о 2 
генеральной дисперсии ос =0.16 и трем измерениям: 
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х= 7.23 х5 =7.8; хз = 7.6;п =3. 
Возьмем уровень значимости а = 0.05. Определим 
7.2+7.8+7.6 
=1.99 





= с = 40/16 =0,4; ира = 1.96, т.к. вероятность 
3 е 
р определена через уровень значимости р=1- 4, и щі} = иад; 
— 2 


2 
иу =и или имеем и, = и, „х . Тогда 
7 то 


2 


х- -И <а<х+ 
а 
1-4% 


< а 
Уп Ја 1%’ 
О оба О ОЕ 


УЗ УЗ 
1.070 =67.99; 
В этом примере приняли случайную величину х с нормальным 


о 
законом распределения с параметрами а и —=. Благодаря этому 
Уп 
а= х 
квантиль и = ——— үп имеет стандартное нормальное распределение и 


с 
с вероятностью р=1-а и удовлетворяет неравенству 


—и <и<и " 
1-92 9 


2 
Однако в расчетах вместо с берут, как правило, выборочную 
2 
дисперсию 5’. Это означает, что вместо квантили и нужно 


а-х 
рассматривать величину { = 5. Уп : 


При больших и> 50 величина 5 мало отличается от с и 
поэтому значения и совпадают с #. 

При малых же объемах выборок различие между ѓи и может 
оказаться существенным и, кроме того, само распределение величины 
і не является нормальным. 

Используя общие законы теории вероятности, можно вывести 
формулы, описывающие распределение величины #1. Это 
распределение называется 1-распределением или распределением 
Стьюдента. Это распределение зависит только от степени свободы #, 


по которым подсчитывается дисперсия 5”. Если объем выборки равен 
п, то Ё =п-1. 

Закон распределения Стьюдента с п-1 степенями свободы 
имеет плотность 


ЛӘ 





где Г- известная гамма - функция. 








== — > 
-5 -4 -3 -2-1012345 і 


Рисунок 5.1 — График плотности # - распределения 


4 
= 


Свойства плотности хорошо видны на графике плотности г — 
распределения (рис. 5.1). Они напоминают плотность нормального 
распределения. При / > оо (п > о) выборочная дисперсия стремится 
к генеральной (5° с?) и распределение Стьюдента сближается с 
нормальным. При малых # распределение г сильно отличается от 


нормального. Поэтому роль г — распределения существенна в 
статистике малых выборок, еще называемой микростатистикой. Будем 
обозначать через ғ, квантили 1 - распределения. 


При доверительной вероятности р=1!1-@ для величины Ѓ 
получаем доверительную оценку —# , 
1- 


<1<1 „. Учитывая, что 
1— 


2 








вет а-х 
і = “п, получим ЕРЕ Мп <! я ИЛИ 
5 1-2 45 Е 
2 2 
А 5 
а ЕЕ одета а ЕРЕ 
ра 7а ит 


Эта оценка очень похожа на ранее рассматриваемые оценки, 
только здесь вместо с берется 5 и поэтому вместо и „ необходимо 


брать ‚ вы 
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Значения ѓ „ для различных чисел степеней свободы } и 
1 


2 
уровней значимости © приводятся в соответствующих книгах по 
математической статистике. 

Распределение Стьюдента позволяет оценить математическое 
ожидание при неизвестной генеральной дисперсии. При этом число 
наблюдений может быть очень малым. 

5.1.2 Оценка генеральной дисперсии 

Известно, что М [52 |= с2. Поэтому, если известно 
распределение выборочной дисперсии 5”, то можно оценить и 
генеральную дисперсию д”. 

Распределение величины 5” можно получить с помощью 
распределения Пирсона, которое еще называют 2 — 
распределением (хи — квадрат распределение). 

Для выборки с элементами х,,х....х„ через у’ обозначается 


сумма 
2 


т(х —х 
=» = (9:10) 
і=1 О’ 


Число степеней свободы здесь также равно Г =и- 1 (одно 
ограничение х). Плотность ж?-распределения зависит только от 
числа степеней свободы } (рис. 5.2). 








24 





0 10 20 30 
Рисунок 5.2 — График ғ? -распределения 
Поскольку х? >20, то и плотность расположена на промежутке 


[0, оо] Е 


Кривые асимметричны, но степень асимметрии уменьшается 
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при увеличении }. 


При уровне значимости © доверительная оценка величины у? 
2 2 2 
имеет вид Хо <Х’<Х д. 
2 =” 
2 2 


Можно усмотреть связь между у” и 56 учетом (5.1) можно 
1 РА 2 2 2 


Е > 
записать 5” = (х, -х) = А откуда 7’ = 


п-1 і=1 


Отсюда с вероятностью р = 1 2 справедливо неравенство 





2 


а 
П 2; 
е 2 Я а Ха 
5 2 
2 Ў 
Квантили у; и х „ даются в таблицах. 


2 2 


5.2 Оценка параметров модели 
Представим модели в следующем виде: 


п=1п(х, В), 
где х — факторы (входные величины); 8 — неизвестные 
параметры (коэффициенты); 7 — реакция системы (выходная 


величина). 


Целью анализа экспериментальных данных является 
определение оценок неизвестных параметров р в некоторой заданной 


области факторного пространства Х . Рассмотрим статистическую 
модель (рис. 5.3). 








Процесс И > 








Рисунок 5.3 — Статистическая модель 


В реальных условиях, из-за наличия помехи &, экспериментатор 
измеряет величину у вместо истинного значения выходной величины 
7. Следовательно, опираясь на результаты измерения, нельзя 
получить абсолютно точных значений д. Вместо истинных 
параметров В приходится использовать случайные величины. 
Обозначим эти величины Ё и назовем их оценками |. Тогда, 
оцениваемое уравнение для модели будет иметь вид: 
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Ү=Ү(х,Ь) (2.2) 


5.2.1 Общие требования, предъявляемые к оценкам параметров 

Чтобы правильно и точно оценить параметры модели, оценки 
должны быть: несмещенными, состоятельными, эффективными и 
достаточными. 

Оценки р являются несмещенными, если их математические 
ожидания равны истинным значениям параметров: 

МІЬЈ= В. 

Это значит, что в процессе вычисления параметров модели не должны 
возникать статистические ошибки. 

Оценка называется состоятельной, если при увеличении числа 
наблюдений п до бесконечности она сходится по вероятности к 
истинному параметру: 


вт Р[(Ь– 8) 2=]=0, при любом &>0. 
по 


Достаточное условие для этого 
іт М[Ф- 8) 1=0 


Оценки будут эффективными, если они позволяют получить 
максимальную информацию из наблюдений. Часто бывает, что из 
исследуемого параметра можно найти несколько состоятельных 
оценок. Чтобы выбрать одну из них сравнивают дисперсии всех 
оценок и по минимуму дисперсии получают оценку, которая и будет 
эффективной 


БЫ< р], 


где [Б] - дисперсия оценки ЁБ, Э[р ] - дисперсия любых других 
несмещенных оценок 6. 

Пример. Дано п наблюдений случайной величины Х. Возникает 
вопрос, какую величину принять за оценку математического 
ожидания: среднее выборки или медиану. Известно, что величина Х 


распределена по нормальному закону с дисперсиями: О[х]= о’ и 


п 
Л 2 =: .. 
0[х,]= 0, где х — среднее выборки; ху — медиана; п — объём 
5 2п их 
2 
о В 
выборки; с“ — дисперсия генеральной совокупности Х. Так как 


0[х]< О[х; ], то оценка средней выборки будет эффективной. 
2 
Критерии несмещенности и эффективности следует 
рассматривать одновременно. Может оказаться, что смещенная 
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оценка с меньшей дисперсией будет более предпочтительной, чем 
несмещенная оценка с большей дисперсией. 
5.2.2 Методы оценивания параметров 

Существует несколько различных методов оценивания 
параметров: 

® максимального правдоподобия; 

• моментов; 

• оценивание по Байесу; 

• наименьших квадратов. 

Метод максимального правдоподобия базируется на 
использовании априорной информации, полученной из эксперимента. 
Получают выборку значений случайной величины Х(х,х›,...,х,). 
Рассматривают оцениваемые параметры 6 как случайные величины с 
некоторым законом распределения вероятности. Затем это 
распределение перестраивается таким образом, чтобы получить 
апостериорное распределение вероятности, плотность которого несет 
информацию о возможных значениях Д на основе экспериментальных 
данных Х. Этот метод приводит к эффективным и состоятельным 
оценкам, однако оценки могут быть смещенными. 

Метод моментов является одним из наиболее старых методов. 
При его использовании вычисляются первые и моментов случайной 
величины, которые затем приравниваются выборочным моментам. 
После этого находят п значений оцениваемых параметров д. 

Оценивание по Байесу как и метод максимального 
правдоподобия основывается на использовании априорной 
информации. Определяется плотность распределения вероятностей х, 
и на основе апостериорной информации принимается решение. 

Метод наименьших квадратов (МНК) является самым 
распространенным методом при оценивании параметров модели. 
Поэтому рассмотрим его более подробно на примере линейной 
модели с одной независимой величиной. 

Уравнение модели с одной независимой величиной имеет вид: 


п = Во + Ах (5.3) 
Оценкой уравнения (5.3) будет: 


Уравнение (5.3) на плоскости представляет теоретическую 
линию регрессии, а (5.4) эмпирическую линию регрессии (рис. 5.4). 
Коэффициенты Ро и Рр; являются оценками истинных коэффициентов 


Вои В. 
На рисунке 5.4 обозначены точки (ух) - одно измерение; 
(7;,х;) — выборочное среднее наблюдение при х; (3;,х) — 


предсказанное значение выходной величины у; при х; (7;,5;) — 
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истинное значение выходной величины 7; при х;. Для несмещенных 
оценок 77; =М [У], т.е. 7; есть математическое ожидание у, прих,. 


По результатам опыта вычисляются коэффициенты Во и р; . Если 
бы все экспериментальные точки оказались на теоретической линии 
регрессии, то ӯ; —77; = 0 или 

у; — ОВ) + Вх, ) = 0, 1=1,2,...п (5.5) 
и тогда коэффициенты д и д; могли быть определены решением 
системы уравнений (5.5). 

Однако, в реальных условиях левая часть (5.5) отличается от 

нуля на величину & 


ж 


у; – (Во + Вух;) = е, (5.6) 


А К (0, (РХ, ) 





Рисунок 5.4 — Линии регрессии однофакторной модели 
а — теоретическая линия регрессии; 6 — эмпирическая линия 
регрессии. 


Величина & называется невязкой. Она может быть вызвана 
ошибкой эксперимента или неправильным выбором линейной модели. 
Поэтому возникает задача найти такие коэффициенты уравнения 
регрессии, при которых невязка будет минимальной. Лучшей оценкой 


п 

является выражение Е, =шш. Это выражение приводит к методу 
і=1 

наименьших квадратов: 


Ф-> Ро -п)' = УРО А-х) (5.7) 


где Р; — число повторных измерений величины у при данном 
значении х;. 

Минимум функции Ф достигается при одновременном 
равенстве нулю частных производных этой функции по всем искомым 
коэффициентам: 
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ОФ ОГУ В - Въ Вих)? _ 
ОВ 980 

әр УВ; - В - Вх)? а Е 
28 08, =-2>.В (5, - Вә – Вх,)х; =0 


После замены бә и б; их оценками Бо и Ё; получаем систему 
нормальных уравнений: 


>В; – Во = Вх; ) =0 





-2> РО, – В – Віх;) =0 
(5.8) 





ИЛИ 
Ро, – Во = Віх)х; = 0 
УВӯ, = У, Вь +), Вых; 59 
(5.9) 
УВӯх = У, Рх +УРЫх 


Решая систему нормальных уравнений относительно Ро и Ё; 
получаем: 





Вә р, г: БАУ – У Вх: У Руд, (5.10) 
УВУ Вх - ОВ)" 
И: ВРУ, (5.11) 





Ув вх - (Вх) 


5.3 Проблема обеспечения точности и достоверности результатов 

Проблема обеспечения точности и достоверности результатов 
можно решить следующим образом. Обозначают точность оценки 
величиной Е. Из-за наличия стохастичности ограниченности числа 


опытов № получают оценку Ё . Задаются неравенством 





Е-Е <= (5.12) 
и вероятностью Р, что неравенство (5.12) выполняется 
О=Р\Е- <) (5.13) 
где О — называют достоверностью оценки. Вводят 
относительную точность оценки 
& 
20 ЛА (5 . 14) 
Е 
с учетом которой достоверная оценка будет равна 
гоа 
о) о (5.15) 








Если известен закон распределения [Е – Е 





‚ то с помощью 


анализа формулы (5.13) или (5.15) можно определить количество 
реализаций №. 
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В тех случаях, когда закон распределения |Е -Ё найти не 


удается, то выдвигают предположение о характере закона 
распределения случайной величины Ё. 

Рассмотрим взаимосвязь точности и достоверности, когда в 
качестве показателей эффективности Е, выступают вероятность Р, 
математическое ожидание а и дисперсия с 

Пусть вероятность появления некоторого события А, которое 
определяется состояниями процесса функционирования исследуемой 
системы, равна 

Р=Р(А) 
В качестве оценки вероятности Р в данном случае выступает 
Б=”. 
№ 

где т - число положительных исходов. Тогда соотношение 
(5.13), связывающее точность и достоверность оценок с количеством 
реализаций, будет иметь вид 


т 
Р-Р – — <= О, 
М = 





РР-е<"<Р+е=0 (5.16) 
М 

Для ответа на вопрос о законе распределения величины Ё 
представим эту частность в виде 

Е 

М №. 
так как количество наступлений события А в данной реализации из № 
реализаций является случайной величиной &, принимающей значения 
х1=1 с вероятностью Р и х›=0 с вероятностью (/-Р). 

Математическое ожидание и дисперсия случайной величины & 
будут равны 


1“ 


м[|=хР+х›(1-Р)=1.Р+0(1-Р)=Р, 

р(]= - м) Р+[ -м (8) Р)= 

=(-Р)?Р+(0-Р)’(-Р)=РА-Р) 
Тогда 


мм - м >» - мм] Р 


Это соотношение говорит о несмещенности оценки Р для 
вероятности Р. 


С учетом независимости значений величины х, получаем 
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Є т 1 м 1 
РР] | р 2,5 -==№РЕ]- 


1 
=— Р-Р) 
М 
В силу центральной предельной теоремы теории вероятностей 
т 
при достаточно больших № можно рассматривать частность № как 


случайную величину, с нормальным законом распределения 
вероятностей с математическим ожиданием Р и дисперсией 


Р(1- Р). 


С учетом квантиля нормального распределения вероятностей 


ір точность оценки Р 


Ег Р-Р) (5.17) 


ео 


Количество реализаций для получения оценки Р с точностью = 
и достоверностью О будет равно 


1 
№=Р Ра-Р).— (5.18) 
0 о 
Е 
1+0 
Квантиль 1, порядка ф= о. находится из специальных 
таблиц. 
Пример. Необходимо рассчитать количество реализаций № при 
статистическом моделировании системы 5 , когда в качестве 


показателя эффективности используется вероятность Р при 
достоверности О = 0.95 (1, =1.96) и точности = = 0.01;0.02;0.05. 


Ввиду того, что значения Р до проведения эксперимента 
неизвестны, то вычисляют множество оценок № для диапазона 
возможных значений Р, т.е. от 0 до 1 с дискретностью 0.1. Результаты 
расчетов с использованием выражения (5.18) представлены в таблице 
5.1. 

Из таблицы 5.1 видно, что при переходе от Р=0.1 (0.9) и 
Р=0.5 количество реализаций № возрастает примерно в 2,5 раза, а 
при переходе от ==0.05и ==0.01 количество реализаций возрастает 
примерно в 25 раз. 


Таблица 5.1 
Результаты расчетов 


Вероятность Р Точность = 
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При тактическом планировании машинных экспериментов, 
когда значение Р неизвестно, поступают следующим образом. Берут 


произвольно значение №, определяют по формуле Р, = > ‚а затем по 
0 

формуле (5.18) проводят вычисления, в которой вместо Р подставляют 

Ро. Такая процедура оценки № может выполняться многократно. 

Если отсутствует возможность получения каких-либо 
априорных сведений о вероятности Р, то в таких случаях 
целесообразно задавать относительную точность результатов 
моделирования є. 


Для этих случаев формула (4.78) принимает вид 
2 | 2 

М№=12 Р-Р) =2? (1-Р). 

Ф 2502; Фф 2 
о ЕЕ 
Оценку показателя эффективности Ё модулируемой системы 
можно дать по результатам определения среднего значения некоторой 
случайной величины. Если случайная величина $ имеет 


2 . 9 
математическое ожидание а и дисперсию о и принимает в гтои 
реализации значение х, то в качестве оценки математического 








(5.19) 


ожидания а принимают 


При больших значениях № в силу предельной центральной 
теоремы теории вероятностей среднее арифметическое х будет иметь 
распределение, близкое к нормальному с математическим ожиданием 
а и дисперсией с /№. 

Точность оценки для математического ожидания равна 


О 
Р Ум 


А ‚ а количество реализаций 


(5.20) 


Если же в качестве показателя эффективности Е выступает 
2 
дисперсия о’, а в качестве ее оценки используется выборочная 
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2 
дисперсия $,то математическое ожидание и дисперсия будут 


м[8°]-(%-1)°: 


д5], 


где ШИ, — центральный момент четвертого порядка случайной 
величины. 
Для частного случая, когда случайная величина имеет 
нормальное распределение, когда и, = 30, получаем 
р 
2 20“ _ 21 


р = 
Е” в, 


№=1 (5.21) 


2 
Для дисперсии о’ точность оценки 





отсюда количество реализаций 
Е 2 Е 4} у 2 о 1 
М№=1, А О Є или № = 1, |5: 
өд 20 
Из выражений (5.19-5.21) видно, что количество реализаций 
существенно зависит от дисперсии оцениваемой случайной величины. 


Поэтому выгодно выбирать те оцениваемые показатели 
эффективности Е, которые имеют малые дисперсии. 


5.4 Проблема уменьшения дисперсии оценок 

Существуют методы, с помощью которых можно увеличить 
точность при заданном числе реализаций и, наоборот, при заданной 
точности оценок сократить число реализаций. Эти методы используют 
априорную информацию о структуре и поведении моделируемой 
системы. 

В качестве такого метода можно использовать метод 
коррелированных реализаций (выборок), который используется в 
задачах сравнения двух или более альтернатив, например, можно 
сравнивать две системы 51 и 52. 

Пусть для систем $; и 52 оценки ёи 2, имеют дисперсии 


Ба | р[2,] и коэффициент корреляции оценок 42,,2, равен К[а,а,|, 


то дисперсию погрешности оценки 4 = @ – 25 разности А = ау -а> 
можно найти из соотношения 
Ба |= оја, —а, |= Ра ]+ 0[2,|-2К[,2,]. оо», (5.22) 
где ОТ = Ба, | О — Ба, | 


При независимом моделировании вариантов системы с 
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использованием различных реализаций псевдослучайных 
последовательностей коэффициент корреляции оценок равен 
К [2 ‚4 ] = 0, 
р,141= [а |+ 51, | (5.23) 
следовательно, Ба 1<р,, [а ],т.к. при имитационном моделировании 
удается получить А > 0. 


5.5 Регрессионный анализ 

После нахождения коэффициентов модели возникает задача 
установить пригодность модели и значимость коэффициентов. С этого 
момента метод наименьших квадратов превращается в регрессионный 
анализ. Применение регрессионного анализа возможно только при 
выполнении следующих предположений. 

1. Математическое ожидание величины у, при заданном 
значении х,является линейной функцией по параметрам, т.е. модели 
должны быть линейными по параметрам. 


2. Значения х; Не ЯВЛЯЮТСЯ случайными величинами. 


а й 2 
3. Дисперсия ошибки ғ, равна дисперсии величины у,: 0; =Оу 


. Эта дисперсия может быть постоянной либо зависимой от х. 

4. Различные измерения величины у взаимно независимы. 

При выполнении этих четырех условий МНК дает несмещенные 
оценки Ро и Б; параметров буи Ду. 

В случае нахождения доверительной области для 
коэффициентов бои б; должно выполняться еще одно предположение. 

5. Условие распределения у; при заданном значении Хх; 


нормально относительно математического ожидания 177; = М[ У, ]. 
5.6 Проверка адекватности модели 
Для того чтобы принять решение относительно модели 
необходимо проверить гипотезу: линейная модель по параметрам 
удовлетворительно описывает экспериментальные данные. 
Проверку адекватности модели выполняют с помощью 
критерия Фишера: 


И 
= г 5 › 
е 
где 52 —– дисперсия неточности модели, 8. —– дисперсия ошибки 


эксперимента. 
Если Е< А „, ТО гипотеза о том, что линейная модель 


адекватна, принимается. Здесь Р_„ —— табличное значение; о 
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уровень значимости; Е — вычисленное значение. 

При Е>Е_„ гипотеза адекватности линейной модели по 
параметрам отвергается. 

Обратимся к рисунку 5.4. Дисперсия 52 характеризует ошибку 
эксперимента, которую можно получить путем возведения в квадрат 
разности (уу – у;)с последующим суммированием по всем 


экспериментальным точкам: 


УО У, 
52 жа (5.24) 


е п 
УР, 
1= 


Знаменатель дисперсии равен разности между общим 
п 
количеством экспериментальных точек УР; и числом наложенных 
1=1 
связей п, т.к. каждая выборка дает одну связь. 
Дисперсия $2 характеризует неточность подгонки модели, её 





неадекватность. Поэтому её можно получить из разности (у – ӯ,), 
которая после возведения в квадрат и суммирования представляет 
собой сумму квадратов отклонений относительно эмпирической 
линии регрессии: 





бы (4.23)(5.25) 


Здесь число степеней свободы равно (п-2). Две связи вызваны 
двумя ограничениями, связанными с определением ри В. 


Если математическое ожидание дисперсии 5” равно дисперсии 
выходной величины 77: 
М[52]=0о?2, 
то полученная эмпирическая модель будет корректной. 
5.7 Определение дисперсии воспроизводимости эксперимента 
Дисперсия 5 " характеризует ошибку эксперимента и называют 
дисперсией воспроизводимости эксперимента. Она может быть 


получена путем объединения дисперсий 5 _ и 5 С 

ОУ ВО. Она 

5, а 1=1 ј=1 - і=1 Е = (5.26) 
(УР, -п) +(п-– 2) Р:=2 
і=1 і=1 
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2 _ с2 
Более простое определение 5, = 5;. 


5.8 Проверка однородности дисперсий 
Простейшим критерием проверки однородности дисперсий 
является критерий Фишера, который представляет из себя отношение 
большей дисперсии к меньшей: 


2 
_ Этах 


52 


тіп 


Е 





Полученная величина Ё сравнивается с табличной величиной 
Ро. Если Ё> Е;_„, то дисперсии значимо отличаются друг от друга и 


тогда они неоднородны. Если сравнивается много дисперсий (более 
двух) и одна из дисперсий значимо превышает остальные, то можно 
применять критерий Кохрена — это отношение максимальной 
дисперсии к сумме всех дисперсий: 


ш 
52 
1=1 


Гипотеза об однородности подтверждается, если С<С,, где С, — 
табличное значение. 


5.9 Принятие решений после построения модели процесса 

Остановимся, в начале, на принятии решения, когда линейная 
модель адекватна. 

При этом возможны три варианта. 

1. Все коэффициенты регрессии значимы. 

2. Часть коэффициентов регрессии значима, а часть незначима. 

3. Все коэффициенты регрессии незначимы. 

Для первого варианта могут быть приняты следующие решения: 
окончания исследования, переход к планам второго порядка и 
движение по градиенту. Переход ко второй серии планирования для 
модели второго порядка позволяет получить экстремум параметра. 

Во втором варианте приходится выдвигать гипотезу, что 
влияние данных эффектов незначимо. Причиной незначимости 
коэффициентов может быть неудачный выбор интервалов 
варьирования и включение второстепенных факторов, большая 
ошибка опыта и т.д. 

Для этого либо расширяют интервалы варьирования по 
факторам, либо увеличивают число параллельных опытов. 

И последний вариант, линейная модель адекватна, все 
коэффициенты не значимы. Это может произойти вследствие большой 
ошибки эксперимента или узких интервалов варьирования. 
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Контрольные вопросы 

1. В чем заключается оценка математического ожидания? 

2. В чем заключается оценка генеральной дисперсии? 

3. В чем заключается оценка параметров модели? 

4. Перечислите общие требования, предъявляемые к оценкам 
параметров. 

5. Какие методы оценивания параметров Вы знаете? 

6. В чем заключается проблема обеспечения точности и 
достоверности результатов? 

7. В чем заключается проблема уменьшения дисперсии 
оценок? 

8. Каким образом проводят проверку адекватности модели? 

9. Как определилить дисперсию воспроизводимости 
эксперимента? 

10. В чем заключается проверка однородности дисперсий ? 

11. Какие решения могут быть приняты после построения 
модели процесса? 
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6 Математические схемы описания технических систем 


6.1 Непрерывно-детерминированные модели (О — схемы) 

Математические схемы такого вида отражают динамику 
процессов, протекающих во времени в системе. Поэтому они 
называются О — схемы. Частным случаем динамических систем 
являются системы автоматического управления. 

Линейная автоматическая система описывается линейным 
дифференциальным уравнением вида 


а"у Я а"х а"! х 
т ОЕ О +В 


где х(1) — задающее воздействие или входная переменная 





ЕРЫ (6.1) 


т т-1 
ГА 


системы; у(1) — состояние системы или выходная переменная; а; ‚В; 


— коэффициенты; / — время. 
На рисунке 6.1 представлена укрупненная функциональная 


схема системы автоматического управления, где Һ (0) = х(ї) – у(ї) — 


сигнал ошибки; р" (Г) — управляющее воздействие; /(1) — 


возмущающее воздействие. Данная система основана на принципе 
отрицательной обратной связи, так как для приведения выходной 
переменной у(1) к ее заданному значению используется информация 
об отклонении между ними. По ней можно разработать структурную 
схему и математическую модель в виде передаточной функции или в 
виде дифференциального уравнения (6.1), в котором для простоты 
предполагается, что точки приложения возмущающих воздействий 
совпадают с входом системы. 





КУ 
х( Управляющая Объект 00 
х п) часть һ"() | управления 

















Датчик 


Рисунок 6.1 — Структура системы автоматического управления 


Непрерывно — детерминированные схемы (0 — схемы) 
выполняются на аналоговых вычислительных машинах (АВМ). 
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6.2 Дискретно-детерминированные модели (Е — схемы ) 

Основным видом дискретно-детерминированных моделей 
является конечный автомат. 

Конечным автоматом называют дискретный преобразователь 
информации, способный под воздействием входных сигналов 
переходить из одного состояния в другое и формировать сигналы на 
выходе. Это автомат с памятью. Для организации памяти в описание 
автомата вводят автоматное время и понятие состояние автомата. 

Понятие «состояние» автомата означает, что выходной сигнал 
автомата зависит не только от входных сигналов в данный момент 
времени, но и учитывает входные сигналы, поступающие ранее. Это 
позволяет устранить время как явную переменную и выразить 
выходные сигналы как функцию состояний и входных сигналов. 

Всякий переход автомата из одного состояния в другое 
возможен не ранее, чем через дискретный интервал времени. Причем 
сам переход считается, происходит мгновенно, то есть не учитывают 
переходные процессы в реальных схемах. 

Существует два способа введения автоматного времени, по 
которому автоматы делятся на синхронные и асинхронные. 

В синхронных автоматах моменты времени, в которых 
фиксируются изменения состояний автомата, задаются специальным 
устройством — генератором синхросигналов. Причем сигналы 
поступают через равные интервалы времени — Д!. Частота тактового 
генератора выбирается такой, чтобы любой элемент автомата успел 
закончить свою работу до появления очередного импульса. 

В асинхронном автомате моменты перехода автомата из одного 
состояния в другое заранее не определены и зависят от конкретных 
событий. В таких автоматах интервал дискретности является 
переменным. 

Также существуют детерминированные и вероятностные 
автоматы. 

В детерминированных автоматах поведение и структура 
автомата в каждый момент времени однозначно определяются 
текущей входной информацией и состоянием автомата. 

В вероятностных автоматах они зависят от случайного выбора. 

Абстрактно конечный автомат можно представить как 
математическую схему (Е — схему), которая характеризуется шестью 
видами переменных и функций: 

1) конечное множество х(1) входных сигналов (входной 
алфавит); 

2) конечное множество у(1) выходных сигналов (выходной 
алфавит); 

3) конечное множество 5(1) внутренних состояний (алфавит 
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состояний); 

4) начальное состояние автомата зо, 0 ЕС; 

5) функция переходов ф(<,х) автомата из одного состояния в 
другое; 

6) функция выходов Ҹ(2, х) автомата. 

Абстрактный конечный автомат имеет один вход и один выход. 
В каждый дискретный момент времени 1=0, [,2,... Е — автомат 
находится в определенном состоянии 5(1) из множества 2 — состояний 
автомата, причем в начальный момент времени /=0 он всегда 
находится в начальном состоянии <(0)=о . В момент 1 , будучи в 
состоянии <(1), автомат способен воспринять на входном канале 
сигнал х()єХи выдать на выходном канале сигнал 
уб) = 107), х(2)1, переходя в состояние 

:10—0+1) = ф[2(7), х(0)] , где 2(0) є 4, х(7) є Х. 

Абстрактный конечный автомат реализует некоторое 
отображение множества слов входного алфавита Х на множество слов 
выходного алфавита У , то есть, если на вход конечного автомата, 
установленного в начальное состояние <о , подавать в некоторой 
последовательности буквы входного алфавита х(0), х(1),2х(2),... , 
которые составляют входное слово, то на выходе автомата будут 
последовательно ПОЯВЛЯТЬСЯ буквы выходного алфавита 
у(0), у(1), у(2),... образуя выходное слово. 

Следовательно, работа конечного автомата происходит по 
следующей схеме: на каждом ї — ом такте на вход автомата, 
находящегося в состоянии 5(1), подается некоторый сигнал х(і), на 
который автомат реагирует переходом на (1+/)~ ом такте в новое 
состояние (1+1) и выдачей некоторого выходного сигнала. 

В зависимости от способа определения выходного сигнала 
синхронные абстрактные конечные автоматы подразделяются на два 
типа: 

— Е — автомат первого рода, также называется автомат Мили: 


5-1) = ф[2(7), х(2)],7 = 01,2,...; 
уб) = Ф[(7), х(0)],7 = 04,2,...; 
— Е — автомат второго рода: 
201+ 1) = ф[2(2), х(7) 1,7 = 01,2,...; 
у) = Ч [< х0 0], = 12,3,...; 


Автомат второго рода, для которого 
уб) = Ч <): = 0/,2,... 
называется автомат Мура — функция выходов не зависит от 
входной переменной х(1). 
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Чтобы задать конечный Е — автомат, необходимо описать все 
элементы множества Ё = (2, х, у,ф, Ф, <0). 


Существует несколько способов задания работы Е — автоматов 
среди которых наибольшее применение нашли табличный, 
графический и матричный. 


6.3 Дискретно — непрерывные модели 
Процессы в линейных импульсных и цифровых системах 
автоматического управления описываются дискретно — разностными 
уравнениями вида: 


ВАК уп) + ВЕ 1А у(п) +...+ Боу(п) = х(п), (6.2) 


где х(и) — решетчатая функция входного сигнала; у(и) — 
решетчатая функция выходного сигнала, которая определяется 


решением уравнения (6.2); р; — постоянные коэффициенты; АК у(п)– 
разность к — го порядка; ѓ=пТ, где иТ – п-ый момент времени, Т — 
период дискретности ( в выражении (6.2) период условно принят за 
единицу). 

Уравнение (6.2) можно представить в другом виде: 

ару(п = К) +ар руп -к+1) +... +аоу(п) = х(п) (6.3) 

Уравнение (6.3) представляет собой рекуррентное соотношение, 
которое позволяет вычислить любой (1+ /)—й член последовательности 
по значениям предыдущих её членов /,/-/,... и значению х(1+/). 

Основным математическим аппаратом моделирования 
цифровых автоматических систем является 7, — преобразование, 
которое базируется на дискретном преобразовании Лапласа. Для этого 
необходимо найти импульсную передаточную функцию системы, 
задаться входной переменной и, варьируя параметрами системы, 
можно найти лучший вариант проектируемой системы. 


6.4 Дискретно — стохастические модели (Р — схемы) 

К дискретно — стохастической модели относится 
вероятностный автомат. В общем, виде вероятностный автомат 
является дискретным потактным преобразователем информации с 
памятью, функционирование которого в каждом такте зависит только 
от состояния памяти в нем и может быть описано статистически. 
Поведение автомата зависит от случайного выбора. 

Применение схем вероятностных автоматов имеет важное 
значение для проектирования дискретных систем, в которых 
проявляется статистически закономерное случайное поведение. 

Для Р — автомата вводится аналогичное математическое 
понятие, как и для Е — автомата. Рассмотрим множество С, 
элементами которого являются всевозможные пары (х; ‚<, ), где х; и =, 
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элементы входного подмножества Х и подмножества состояний 7 
соответственно. Если существуют две такие функции Ф и Ч, что с их 
помощью осуществляется отображение С > 7 и С —> У, то говорят, 
что = (4, Х,Ү,Фф, Р) определяет автомат детерминированного типа. 


Функция переходов вероятностного автомата определяет не 
одно конкретное состояние, а распределение вероятностей на 
множестве состояний (автомат со случайными переходами). Функция 
выходов также есть распределение вероятностей на множестве 
выходных сигналов (автомат со случайными выходами). 

Для описания вероятностного автомата введем в рассмотрение 
более общую математическую схему. Пусть Ф — множество 
всевозможных пар вида (3, ,у;), где у, — элемент выходного 
подмножества У. Далее потребуем, чтобы любой элемент множества С 
индуцировал на множестве Ф некоторый закон распределения 
следующего вида: 

элементы из Ф ... (21, У1 (51, У2)...(< к, уу-1)(к› Уу) 

а) ЕН Ру Ру 6 Бе 3-1) Бру 

где рр; — вероятности перехода автомата в состояние гк и 
появления на выходе сигнала у; , если он был в состоянии г; и на его 
вход в этот момент времени поступал сигнал х;. 

Число таких распределений, представленных в виде таблиц 


равно числу элементов множества С. Если обозначить это множество 
таблиц через В, то тогда четверку элементов Р= (7, Х,Ү, В) называют 


вероятностным автоматом ( Р — автоматом). При этом 
кл 
Э; Хь} =; 
К=1 ј=1 
Частным случаем Р — автомата, задаваемого как 


Р=(7, Х,Ү, В) являются автоматы, у которых либо переход в новое 
состояние, либо выходной сигнал определяются детерминировано(й 
—  детерминированный вероятностный автомат У — 
детерминированный вероятностный автомат соответственно). 

Очевидно, что с точки зрения математического аппарата задание 
У — детерминированного Р — автомата эквивалентно заданию 
некоторой марковской цепи с конечным множеством состояний. В 
связи с этим аппарат марковских цепей является основным при 
использовании Р— схем для аналитических расчетов. Подобные Р — 
автоматы используют генераторы марковских последовательностей 
при построении процессов функционирования систем или 
воздействий внешней среды. 

Марковские последовательности, согласно теореме Маркова 
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[9], — это последовательность случайных величин, для которой 
справедливо выражение 
1 М 
іт (— 2. =0, 
Мэ № і=1 
где М — количество независимых испытаний; О — дисперсия. 
Такие Р — автоматы (Р — схемы) могут быть использованы для 
оценки различных характеристик исследуемых систем как для 
аналитических моделей, так и для имитационных моделей с 
использованием методов статистического моделирования. 
У — детерминированный Р — автомат можно задать двумя 
таблицами: переходов (табл.6.1) и выходов (табл.6.2). 


Таблица 6.1 


Таблица переходов 





Таблица 6.2 








Где Р; — вероятность перехода Р — автомата из состояния 2; в 
К 
состояние 2}, при этом >, Ру =1. 
Ј=1 


Таблицу 6.1 можно представить в виде квадратной матрицы 
размерности Кх. Такую таблицу будем называть матрицей 
переходных вероятностей или просто матрицей переходов Р — 
автомата, которую можно представить в компактной форме: 








Ру Во ... Ру 
о в. 
„= 
Ра Ро ... Рек 
Для описания У — детерминированного Р — автомата 


необходимо задать начальное распределение вероятностей вида: 
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И 71 72 256 2-1 7к 


р: Ч 4 5: Чк-1 Чк 
где 4, вероятность того, что в начале работы Р — автомат 
К 
находится в состоянии 2%, при этом »Ах =1. 
К=1 


И так, до начала работы Р — автомат находится в состоянии 740 и 
в начальный (нулевой) такт времени меняет состояние в соответствии 
с распределением О. После этого смена состояний автомата 
определяется матрицей переходов Р. С учетом 27 размерность 
матрицы Р, следует увеличить до (& +1) х(& +1), при этом первая 
строка матрицы будет (4,4, 4ь,...,Ак), а первый столбец будет нулевым. 
Пример. У— детерминированный Р— автомат задан матрицей 
переходов: 
0,50 0 0 0.50 


0 

о о о 100 0 
Р ={0 0 0,75 0 0,25 

о о 0,40 0 0о,60 

0 100 0 0 0 





и таблицей выходов 
Таблица 6.3 
Таблица выходов 


С учетом матрицы переходов граф переходов вероятностного 
автомата представлен на рисунке 6.2. 

Требуется оценить суммарные финальные вероятности 
пребывания этого автомата в состоянии 72 И 73, т.е. когда на выходе 
автомата появляются единицы. 








Рисунок 6.2. — Граф переходов 
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При аналитическом подходе можно использовать известные 
соотношения из теории марковских цепей и получить систему 
уравнений для определения финальных вероятностей. Причем 
начальное состояние можно не учитывать в виду того, что начальное 
распределение не оказывает влияние на значения финальных 
вероятностей. Тогда матрица переходов без первого столбца и строки 
примет вид: 


о 0 100 0 
0 0,75 0 0,25 
0 0,40 0 л 0,60 

1000 о о 


С | 
| 
С: 








где С=(сь)=(с1],со,сз,сд) — финальная вероятность 
пребывания У — детерминированного Р — автомата в состоянии 1. 
В результате получаем систему уравнений: 


ор 64 
Сә = 0,75с, +0,40с; 
6.4 
Ке (6.4) 


сд =0,25с> + 0,60с; 
К данной системе следует добавить условие нормировки: 
С + Со + С3 + Сд =1 (6.5) 
Теперь решая систему уравнений (6.4) совместно с (6.5), 
получаем: 


= 80 (7 ть 
23 23" 23 
Таким образом, при бесконечной работе заданного автомата на 
его выходе будет формироваться двоичная последовательность с 


С 13 
вероятностью появления единицы, равной: со + сз = — = 0,5652. 


Кроме аналитических моделей в виде Р — схем можно 
применять и имитационные модели, реализуемые, например, методом 
статистического моделирования. 


6.5 Непрерывно-стохастические модели (О — схемы) 

Такие модели рассмотрим на примере использования в качестве 
типовых математических схем систем массового обслуживания, 
которые называют О — схемами. Такие О — схемы применяются при 
формализации процессов функционирования систем, которые по 
своей сути являются процессами обслуживания [5, 9,10]. 
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К процессам обслуживания можно отнести: потоки поставок 
продукции некоторому предприятию, потоки деталей и 
комплектующих изделий на сборочном конвейере цеха, заявки на 
обработку информации ЭВМ от удаленных терминалов сети ЭВМ. 
Характерным признаком для функционирования таких систем или 
сетей является случайное появление заявок на обслуживание. Причем 
в любом элементарном акте обслуживания можно выделить две 
основные составляющие: ожидание обслуживания и , собственно, сам 
процесс обслуживания заявки. Представим это в виде некоторого 1-го 
прибора обслуживания П; (рис.6.3), состоящего из накопителя заявок 


Н; , в котором может находится одновременно 4, заявок; К; — канал 


обслуживания заявок. 
На каждый элемент прибора П; поступают потоки событий, в 
накопитель Н; поток заявок 4, , на канал К; — поток обслуживания И;. 


= 











5 


Рисунок 6.3 — Прибор обслуживания 


Потоки событий могут быть однородными, если он 
характеризуется только последовательностью поступления этих 
событий ( О0</ < <...<1,), или неоднородными, если он 


характеризуется набором признаков события, например таким 
набором признаков: источник заявок, наличие приоритета, 
возможность обслуживания тем или иным типом канала и т.п. 

Обычно при моделировании различных систем применительно к 
каналу К; можно считать, что поток заявок 4; на входе К; образует 
подмножество неуправляемых переменных, а поток обслуживания И; 
— образует подмножество управляемых переменных. 

Те заявки, которые по различным причинам не обслуживаются 
каналом К;, образуют выходной поток У;. 

Приборы могут подключаться последовательно и (или) 
параллельно, что предоставляет возможность моделировать процессы 
массового обслуживания для разных объектов. Такой принцип 
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называется агрегатированием. 


6.6 Стохастические минимаксные модели 

Эти модели можно отнести к оптимальным стохастическим 
моделям. 

Во многих случаях при построении модели не все условия 
заранее известны. Эффективность нахождения модели здесь будет 
зависеть от трех факторов: 

® заданных условий ху, Х»...,Хи; 

® неизвестных условий у, у,...,уб 

• зависящих от нас факторов иџ,и›,...,и,, которые необходимо 
найти. 

Показатель эффективности решения такой задачи имеет вид: 


Й = М(х], хо,...,х,; у, У2,:. Ук Щ,и2,..5И) 

Наличие неизвестных условий у; переводит задачу оптимизации 
в задачу о выборе решения в условиях неопределенности. Задача 
становится чрезвычайно сложной. 

Особенно задача осложняется для случаев, когда величины у; не 
обладают статистической устойчивостью, то есть неизвестные 
условия у; нельзя изучить с помощью статистических методов. Их 
законы распределения либо не могут быть получены, либо вовсе не 
существуют. 

В этих случаях рассматриваются комбинации всевозможных 


М В Е ж 
значений Ү:о = (уу, у2,...,ук) таким образом, чтобы получить как 


«наилучшее», так и «наихудшее» сочетания значений переменных у;. 
Тогда в качестве критерия оптимизации рассматривается: 


Я = М(х, с,и) 
х = х(х1,х2,....^и),а = а(@1, 42 ,...@ р) 


и = и(иу,иә,...и) 
Критерий Вальда 
При использовании этого критерия всегда ориентируются на то, 
что случайные факторы сложатся в наихудшую комбинацию. Такой 
критерий называют пессимистическим: 
тах тіп И’(х, си). 


Пример. Имеется четыре стратегии игрока: иъизиз, ид; три 
варианта условий: у, у,уз и известные значения х; 1=1,2,3,4, ]=1,2,3. 

Требуется найти оптимальное решение с использованием 
критерия Вальда. 

Составим матрицу игры наихудших результатов (табл. 6.4). 
Согласно критерия Вальда: 0; = аз = 0,25. 
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Таблица 6.4 


Матрица игры 
о ИМЕН 


Следовательно оптимальной будет стратегия из. 
Критерий Гурвица 
При этом критерии предполагается, что с вероятностью Дд 
сложатся наилучшая, а с вероятностью (1-5) наихудшая комбинации 
случайных факторов: 
тах[ 5 тах УИ (х, <, и) + (1- В)тшИ(х,а,и)] 
и [04 а 





(6.6) 


0<6д<1 

При Д=0 критерий Гурвица переходит в критерий Вальда, при 
В=1 получается так называемый критерий «здорового оптимизма». 
Коэффициент 5 выбирается из субъективных соображений. Если 
ситуация очень опасная, то его выбирают ближе к 1. 

Пример. Берется тот же пример (табл.6.4). Принимается Д=0,4. 
Достраивается таблица 6.4 столбцом, соответствующим оптимальной 
оценке (табл. 6.5). Максимальное значение 773=0,47 соответствует 
стратегии из. 


Таблица 6.5. 
Матрица игры 


у, ИИ 


095-000-952 ооо оа 





По критерию (6.6) получаем: 7; = 0,40; + 0,60;. 


Контрольные вопросы 
1. Для каких объектов применяются непрерывно — 
детерминированные модели (0 — схемы)? Перечислите особенности 
этой схемы. 
2. Для каких объектов применяются дискретно — 
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детерминированные модели (Е — схемы)? Перечислите особенности 
этой схемы. 

3. Для каких объектов применяются дискретно — 
непрерывные модели? Перечислите особенности этой схемы. 

4. Для каких объектов применяются дискретно — 
стохастические модели (Р — схемы)? Перечислите особенности этой 
схемы. 

5. Для каких объектов применяются непрерывно — 
стохастические модели (О — схемы)? Перечислите особенности этой 
схемы. 

6. Для каких объектов применяются стохастические 
минимаксные модели? Перечислите особенности этой схемы. 
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